NAVAL  POSTGRADUATE  SCHOOL 
Monterey ,  California 


w. 


DTIC 


-LECTE 
AUG  1  6  1988 


n 


A  SPARES  OPTIMIZATION 
U.S.  MARINE 

MODEL  FOR  DEPLOYABLE 
CORPS  UNITS 

by 

Paul  R. 

Yor  io 

March 

1988 

Thesis  Advisor : 

Donald  P. 

Gaver 

Approved  for  public  release:  distribution  is  unlimited 


w 

;1 


I 

'i*«| 


J| 


ft 


_  f 


CURITV  CLASS. f.CATiON  O 


i»  report  security  classification 

Unclassified 


2a.  security  classification  authority 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION i  OOWNuRADiNG  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NuMBER(S) 


3  DISTRIBUTION  I  availability  OF  REPORT 

Approved  for  public  release; 
Distribution  is  unlimited 


}  MON,  ORING  ORGANIZATION  REPORT  NUMBERS/ 


6a  NAME  OF  PERFORMING  ORGANIZATION  6o  OFFiCE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

(if  applicable) 

Naval  Postgraduate  School  ,,  ,  ,  0rtCr„ra,alia<.a  q„H 


6c  ADDRESS  ;C/fy,  Stare,  ana  ZIP  Code) 


Monterey,  California  93943-5000 


Naval  Postaraduate  School 


7o  ADDRESS  (Ciry.  State,  ana  ZIP  Code) 

Monterey,  California  93943-5000 


8a  NAME  OP  FuND'NG  SPONSORING 
ORGANIZATION 


8c  ADDRESS  (C/ty.  Stare  ana  ZiPCoael 


1 1  title  (include  Security  Gass/bcenon) 


8o  OFFICE  SYMBOL  9  PROCUREMENT  .NSTRUMENT  IDENTIFICATION  NUMBER 
(if  applicable) 


PROGRAM 
FlEMENT  NO 


UND'NG  NUMBERS 

RROjECT 

NO 

TASK 

NO 

•VOPK  jN'iT 
•:rcss,ON  no 

A  SPARES  OPTIMIZATION  MODEL  FOR  DEPLOYABLE  U.S.  MARINE  CORPS  UNITS 


12  PERSONAL  AuTHOR(S) 


13a  TYPE  OF  REPORT 

Master  1  s  Thesis 


YORIO ,  Paul  R. 


•3d  -  ME  COVERED 
PROM  TO 


14  DATE  OF  REPORT  < Year.  Month,  Day)  IS  PAGE  C',uNT 

1988  March 


16  supplementary  notation  The  views  expressed  in  this  thesis  are  those  of  the  author 
and  do  not  reflect  the  official  policy  or  position  of  the  Department  of 

”  "  Yip n  h  .  _  _ 


18  SuBiECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Single  period  inventory  model,  GAMS, 
spares  optimization 


17 

COSAT!  CODES  1 

FiELO 

GROUP 

SUBGROUP  ] 

’9  ABSTRACT  (Continue  on  reverse  it  necessary  and  identify  by  block  number) 

The  U.S.  Marine  Corps  deploys  Marine  Air-Ground  Task  Forces  ( MAGTFs )  by 
airlift  or  sealift  to  participate  in  numerous  short-term  exercises.  These 
exercises  are  of  such  duration  that  resupply  of  the  MAGTF  by  strategic 
airlift  or  sealift  is  not  practical.  Thus,  only  stocked  spare  parts  are 
available  for  repairs  during  the  exercise. 

A  model  is  developed  which  provides  the  operational  commander  with  a 
stockage  policy  for  spare  secondary  reparables  (e.g.,  tank  engines,  amtrack 
transmissions,  etc.)  that  optimizes  the  probability  of  successful  mission 
completion  subject  to  weight  or  volume  constraints  imposed  by  the  MAGTF ' s 
mode  of  deployment.  Optimization  of  this  stockage  policy  is  stochastically 
modeled  using  data  from  the  Marine  Corps  Integrated  Maintenance  Management 
System  data  base  and  then  solved  as  an  integer  program.  (Continued) 


20  distribution  1  availability  of  abstract 

Q UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Donald  P.  Gaver 


21  abstract  security  CLASSIFICATION 
□  otic  USERS  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  c2 c  OFFICE  SYMBOL 

408-646-2605  55Gv 


DO  FORM  1473,  84  mar 


83  APR  edition  may  be  used  until  enhausted 
aii  other  editions  are  obsolete 

i 


SECURITY  CLASSIFICATION  of  this  PAGE 

•  U.S.  Prtntlnf  0">c«  ISM— 


Block  19  Abstract.  (Continued) 

The  integer  program  is  coded  using  the  Generalized  Algebraic 
Modeling  System  language  and  solved  using  the  Zero/One  Optimization  Methods 
mixed  integer  program  solver.  Operational  data  for  a  Marine  Amphibious 
Unit  yields  an  integer  program  with  190  binary  variables  and  26  constraints 
A  solution  within  0.07%  of  optimality  is  obtained  on  an  IBM  3033AP  computer 
in  3.9  seconds  and  on  a  Zenith  Z-248  personal  computer  in  176  seconds. 


^  ACZi-Jhii. :  ■  r  O' 

!  n  ns  CRA-V 
!  DT'iC  TAB 

|  i J T K J ' ' f  t I  'DCO'j 

.  •  '0  j  tjO;  i 

1 - - 


Avuilabii'ty  Codes 

f  AvdH  and i or 
^•ecial 


I 

i 


Approved  for  Public  Release;  distribution  is  unlimited. 


A  Spares  Optimization  Model  for  Deployable 
U.S.  Marine  Corps  Units 


by 


Paul  R.  Yorio 

Captain,  United  States  Marine  Corps 
BS,  U.S.  Naval  Academy,  1982 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  1988 


Author: 


Approved  by: 


Paul  <R.  Yorio 


Donald  P.  Gaver,  Thesis  Advisor 


ABSTRACT 


The  U.S.  Marine  Corps  deploys  Marine  Air-Ground  Task 
Forces  ( MAGTFs )  by  airlift  or  sealift  to  participate  in 
numerous  short-term  exercises.  These  exercises  are  of  such 
duration  that  resupply  of  the  MAGTF  by  strategic  airlift  or 
sealift  is  not  practical.  Thus,  only  stocked  spare  parts 
are  available  for  repairs  during  the  exercise. 

A  model  is  developed  which  provides  the  operational 
commander  with  a  stockage  policy  for  spare  secondary 
reparables  (e.g.,  tank  engines,  amtrack  transmissions,  etc.) 
that  optimizes  the  probability  of  successful  mission 
completion  subject  to  weight  or  volume  constraints  imposed 
by  the  MAGTF' s  mode  of  deployment.  Optimization  of  this 
stockage  policy  is  stochastically  modeled  using  data  from 

the  Marine  Corps  Integrated  Maintenance  Management  System 

V  Cc>r.jpvtc>- 

data  base  and  then  solved  as  an  integerjj^program. 

The  integer  program  is  coded  using  the  Generalized 

-f (3h\-\S) _ j 

Algebraic  Modeling  System ^language  and  solved  using  the 
Zero/One  Optimization  Methods  mixed  integer  program  solver. 
Operational  data  for  a  Marine  Amphibious  Unit  yields  an 
integer  program  with  190  binary  variables  and  26 
constraints.  A  solution  within  0.07%  of  optimality  is 
obtained  on  an  IBM  3033AP  computer  in  3.9  seconds  and  on  a 
Zenith  Z-248  personal  computer  in  176  seconds.  " 


TABLE  OF  CONTENTS 


I.  INTRODUCTION .  1 

A.  PROBLEM  STATEMENT .  2 

B.  INVENTORY  MODELS  AND  SOLUTION  APPROACHES....  7 

C.  THESIS  OUTLINE . 15 

II.  DATA  ISSUES  AND  ANALYSIS . 17 

A.  DESCRIPTION  OF  DATA . 17 

B.  DATA  ANALYSIS . 22 

III.  MODEL  FORMULATION  AND  DESCRIPTION . 28 

A.  MODEL  DEVELOPMENT . 28 

B.  TRANSFORMATION  TO  A  BINARY  LINEAR  PROGRAM... 40 

IV.  SOLUTION  METHODOLOGY  AND  IMPLEMENTATION . 45 

A.  SOLUTION  APPROACHES  AND  METHODOLOGY . 45 

B .  IMPLEMENTATION . 50 

V.  RESULTS . 53 

A.  MODEL  TESTING . 53 

B.  MODEL  VERIFICATION . 56 

C.  BAYESIAN  AND  EMPIRICAL  BAYES  ANALYSIS . 60 

D.  COMPUTATION  EXPERIENCE . 63 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS . 67 

APPENDIX  A.  LISTING  OF  SOURCE  CODE  FOR  LANGRANGIAN 

RELAXATION . 70 

APPENDIX  d.  SAMPLE  INPUT  FILE . 75 

APPENDIX  C.  SAMPLE  OUTPUT  FILE . 78 


V 


APPENDIX  D.  GAMS  SOURCE  CODE 


t 

< 

* 


< 

< 


v 


•> 


* 

{ 

\ 


< 

a 

u 

u 

D 

a 

o 

e 

c. 

Ui 

C 


! 


5 


I 

? 


* 

* 

J 

j 


I 

C 

K 


LIST  OF  REFERENCES . 

INITIAL  DISTRIBUTION  LIST 


* 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  thanks  and  gratitude  to 
the  following  individuals  without  whom  this  project  could 
not  have  been  completed. 

-  My  wife,  Juzenilde,  for  her  patience  and  continuous 
support . 

-  My  thesis  advisors.  Professor  Donald  P.  Gaver  and 
Professor  Kevin  Wood,  whose  guidance  and  positive  direction 
kept  this  study  on  course . 

-  Professor  Richard  E.  Rosenthal  who  introduced  me  to  the 
Generalized  Algebraic  Modeling  System. 

-  Major  Robert  D.  Larson,  USMC  of  the  SMU  Operations 
Section,  2d  Supply  Battalion,  2d  Force  Service  Support 
Group,  FMF  who  provided  both  the  impetus  for  this  study  and 
the  data. 

-  Alexander  Meeraus  and  Anthony  Brooke  of  the  GAMS 
Development  Corporation,  Washington,  D.  C.  for  providing  a 
courtesy  copy  of  the  Generalized  Algebraic  Modeling  System 
to  the  Naval  Postgraduate  School. 


REPRODUCED  AT  GO V E R N M r  I )  I 


I. 


INTRODUCTION 


The  combat  structure  of  the  United  States  Marine  Corps 
(USMC)  centers  around  the  concept  of  the  Marine  Air-Ground 
Task  Force  (MAGTF)  which  is  a  combined  arms  force  consisting 
of  ground,  air,  and  combat  service  support  forces.  The  size 
of  a  MAGTF  and  its  specific  composition  are  dependent  upon 
the  particular  mission  assigned  and  the  capability  of  the 
opposing  forces.  However,  all  MAGTFs  are  organized  to  take 
maximum  advantage  of  the  combat  potential  inherent  in  a 
rapidly  deployable,  closely  integrated  air-ground  team  under 
the  control  of  a  single  commander.  A  MAGTF  consists  of  four 
basis  elements:  command  element,  ground  combat  element, 
aviation  combat  element  and  combat  service  support  element. 
It  is  the  mission  of  the  combat  service  support  element  to 
provide  all  logistical  support  for  the  ground  and  aviation 
combat  elements.  [Ref.  l:pp.l-4]. 

In  this  study,  a  model  is  developed  which  provides  the 
operational  commander  with  a  stockage  policy  for  spare 
secondary  reparables  (e.g.,  tank  engines,  amtrack 
transmissions,  etc.)  that  optimizes  the  probability  of 
successful  mission  completion  within  the  constraints  imposed 
by  the  MAGTF ' s  mode  of  deployment  (i.e.,  airlift  or 
sealift).  The  problem  of  determining  this  optimal  stockage 
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policy  is  stochastically  modeled  using  operational  data 
resident  in  the  Marine  Corps  Integrated  Maintenance 
Management  System  (MIMMS)  data  base  and  then  solved  as  a 
mathematical  program.  The  remainder  of  this  chapter  will 
include  a  description  of  the  scenarios  which  are  of 
particular  interest  to  the  U.S.  Marine  Corps,  an  explanation 
of  terms,  and  a  definition  of  the  problem.  Additionally, 
justification  for  the  approach  taken  here  as  compared  to 
those  in  standard  inventory  models  is  provided. 


A.  PROBLEM  STATEMENT 

The  U.S.  Marine  Corps  currently  deploys  MAGTFs  by 
airlift  or  sealift  to  participate  in  numerous  short-term 
exercises  such  as  BRIGHT  STAR  in  Egypt,  GALLANT  EAGLE  in  the 
southwestern  United  States,  and  COLD  WINTER  in  Norway. 
These  MAGTFs  range  in  size  from  2,000  man  Marine  Amphibious 
Units  (MAUs)  to  16,000  man  Marine  Amphibious  Brigades 
(MABs),  depending  upon  the  particular  exercise.  These 
exercises  are  of  such  duration  that  resupply  of  the  MAGTF's 
ground  combat  element  by  strategic  airlift  or  sealift  from 
logistics  bases  in  the  continental  United  States  is  not 
practical.  Furthermore,  by  doctrine  [Ref.  l:pp.  1-4],  a 
MAGTF  is  task  organized  to  maximize  combat  power;  therefore, 
the  number  of  mechanics,  technicians,  and  other  logistics 
personnel  together  with  their  test  and  repair  equipment  are 
kept  to  a  minimum.  Thus,  the  quantities  of  spare  repair 


parts  stocked  by  the  combat  service  support  elements  are  all 
that  are  available  to  support  the  requirements  of  the  ground 
combat  element  for  the  duration  of  the  exercise. 

In  addition  to  these  operational  exercises,  the  Marine 
Corps  deploys  Marine  Amphibious  Units  (MAUs)  to  the 
Mediterranean  Sea  as  part  of  the  U.S.  Sixth  Fleet  on  a 
continuous,  six  month  rotating  basis.  These  Mediterranean 
MAUs  remain  aboard  ship  except  when  conducting  exercises. 
Therefore,  the  inventory  of  spares  maintained  by  their 
combat  service  support  element,  called  a  MAU  Service  Support 
Group  ( MSSG )  ,  is  limited  by  the  cargo  capacity  of  the 
assigned  amphibious  ships.  Even  though  resupply  from  the 
continental  United  States  by  air  is  possible,  it  is  desired 
to  reduce  dependency  on  this  costly  method  [Ref.  2:p.  1]  . 
Thus,  similar  to  the  above  exercise  scenarios,  the  MSSG  must 
stock  sufficient  spares,  within  the  capacity  constraints  of 
the  assigned  shipping,  so  as  to  satisfy  the  anticipated 
demands  generated  by  the  ground  combat  element's 
participation  in  exercises  during  the  course  of  the  six 
month  deployment. 

Presently,  a  problem  of  even  greater  concern  to  Marine 
Corps  combat  service  support  organizations  than  either  of 
those  previously  mentioned  is  the  Marine  Corps's  role  in  the 
rapidly  evolving  maritime  strategy  of  the  U.S.  Navy.  As  the 


amphibious  power  projection  arm  of  the  Navy,  the  Marine 
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Corps,  in  conjunction  with  the  Navy,  has  developed  and 
recently  implemented  the  concept  known  as  the  Maritime 
Prepositioning  Force  (MPF).  According  to  doctrine  [Ref.  3], 
a  MPF  consists  of  a  squadron  of  four  to  five  ships  loaded 
with  combat  equipment  and  supplies.  One  such  MPF  squadron 
is  presently  located  in  each  of  the  following  three  ocean 
areas:  Eastern  Atlantic,  Indian  Ocean, and  Western  Pacific. 
When  a  MPF  operation  is  ordered,  a  fly-in  echelon  consisting 
of  16,500  Marines  and  sailors,  comprising  a  Marine 
Amphibious  Brigade  ( MAB )  and  a  Navy  Support  Element,  is 
airlifted  to  a  benign  airfield  in  the  vicinity  of  the 
objective  area  for  link-up  with  equipment  and  supplies 
aboard  the  MPF  squadron  ships.  The  Brigade  can  be  combat- 
capable  and  ready  to  move  to  an  objective  in  five  days  or 
less  and  can  operate  for  thirty  days,  independent  of  any 
strategic  airlift  or  sealift  from  the  United  States,  using 
the  supplies  and  spare  parts  inventories  on  the  ships  [Ref. 
4:p.  11].  It  is  this  last  requirement,  to  operate 
independently  for  thirty  days  using  organic  supplies  and 
inventories,  that  poses  the  greatest  challenge  to  Marine 
Corps  combat  service  support  planners.  However,  as  with 
supporting  short-term  exercises  and  MAU  deployments,  the 
problem  of  supporting  the  MPF  MAB  simplifies  to  one  of 
stocking  adequate  spare  parts  inventories,  subject  to  the 
capacity  constraints  of  the  MPF  squadron  ships,  to  meet  the 
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demand  during  the  MAB ' s  projected  missions.  The  fact  that 
such  demand  may  vary  unpredictably  must  be  confronted. 

At  the  present  time,  to  solve  spare  stockage  problems 
similar  to  those  described  above.  Marine  Corps  combat 
service  support  planners  determine  the  number  of  secondary 
items  (referred  to  in  Marine  Corps'  terminology  as  secondary 
reparables )  to  stock  for  each  principal  end  item  in  the 
particular  MAGTF's  Table  of  Equipment  (T/E)  by  using  either 
last-period  demand  method  or  a  simple  moving  average  of  past 
or  historical  demands  generated  during  similar  exercises  or 
Mediterranean  deployments.  In  the  situation  where  no 
historical  demands  are  available,  such  as  with  the  new  MPF 
MABs,  stockage  levels  are  determined  using  estimated 
replacement  rates  in  combination  with  the  experience  of 
maintenance  and  supply  personnel.  These  techniques  have 
c  proved  inadequate  [Ref.  2:p.  1]  because  they  fail  to 

consider  such  factors  as  changes  to  a  MAGTF's  Table  of 
Equipment  and  variations  in  the  operational  schedule,  both 
of  which  may  occur  from  one  exercise  or  deployment  to  the 
next.  More  importantly,  they  neglect  to  consider  that  the 
demands  placed  on  the  supply  system  are  the  direct  result  of 
equipment  failures. 

To  overcome  these  inadequacies,  it  is  the  objective  of 
this  study  to  develop  a  model  which  provides  the  operational 
commander  at  the  MAU  or  MAB  level  with  a  stockage  policy  for 
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spare  secondary  reparables  that  optimizes  the  probability  of 
successful  mission  completion  within  the  constraints  imposed 
by  the  mode  of  deployment.  The  proposed  model  is 
sufficiently  general  so  as  to  be  adaptable  to  any  of  the 
above  scenarios  and  it  can  accommodate  any  additions  or 
deletions  to  the  MAGTF ' s  Table  of  Equipment  or  variations 
in  mission  duration  among  any  combination  of  principal  end 
items.  Additionally,  the  model  incorporates  a  measure  of 
effectiveness  which  provides  the  operational  commander  with 
an  indicator  of  the  overall  impact  of  the  stockage  policy  on 
successful  mission  completion.  Furthermore,  the  model 
utilizes  the  operational  maintenance  data  resident  in  the 
Marine  Corps  Integrated  Maintenance  Management  System 
(MIMMS)  to  determine  operationally  based  failure  rates. 
Thus,  by  use  of  this  model  a  nearly  optimal  stockage  policy 
for  secondary  reparables  at  the  retail  level  can  be 
determined  for  a  wide  variety  of  inventory  problems 
currently  facing  Marine  Corps  combat  service  support 
planners.  This  will  result  in  a  more  effective  utilization 
of  the  limited  space  and  weight  resources  available  for 
spares'  inventories  on  shipping  or  aircraft  assigned  to 
deploying  Marine  Corps  units. 


B.  INVENTORY  MODELS  AND  SOLUTION  APPROACHES 


The  stockage  problems  confronting  Marine  Corps  planners 
can  all  be  viewed  as  variations  of  the  single  period, 
constrained,  multi-item  inventory  problem.  The  simplest 
such  single-period  model  is  the  classic  newsboy  problem 
[Ref.  5],  alternatively  referred  to  as  the  Christmas  tree 
problem.  The  characteristic  distinguishing  the  models  of 
this  thesis  from  other  standard  inventory  models,  such  as 
the  classical  economic-order-quantity  or  economic-order- 
interval  models  [Ref.  6:pp.  159-289]  is  that  only  a  single 
time  period,  usually  of  finite  length,  is  relevant; 
therefore,  only  a  single  initial  procurement  or  order  is 
made  [Ref.  6:p.  297].  The  previously  described  scenarios 
all  fall  within  the  context  of  single-period  models  because 
the  demand  for  secondary  reparables  is  the  direct  result  of 
equipment  failures  which  occur  at  infrequent  intervals 
during  the  course  of  the  exercise.  Resupply  either  does  not 
occur  or  is  kept  to  a  minimum,  and  thus  it  can  be  ignored. 
Furthermore,  it  is  Marine  Corps  policy  to  permit  only 
functional  unit  replacement  at  the  field  or  organizational 
levels  of  maintenance.  Actual  repair  of  secondary  reparables 
at  the  intermediate  maintenance  level  is  very  limited  and 
for  all  practical  purposes  performed  only  at  the  depot 
level.  Therefore,  since  the  MAGTFs  of  concern  in  these 
scenarios  are  usually  only  authorized  up  to  the  intermediate 


level  of  maintenance,  the  actual  repair  and  return  of 
secondary  reparables  to  these  retail  levels  of  inventory  is 
insignificant  and  can  be  ignored.  This  is  particularly  true 
when  one  considers  that  some  of  the  exercise  durations  are 
of  the  same  order  of  magnitude  as  the  average  repair  time 
for  some  of  the  secondary  reparables.  Thus,  each  scenario 
consists  of  an  exercise  (or  exercises)  which  is  (are)  of 
finite  duration  and  during  which  repair  and  resupply  are 
essentially  nonexistent,  so  the  combat  service  support 
planner  must  make  a  one-time  decision  of  what  to  stock,  not 
how  much  of  each  item  to  procure.  He  must  do  so  within  the 
constraints  of  the  mode  of  deployment.  This  is  analogous  to 
the  single  period,  constrained,  multi-item  inventory 
problem. 

There  are  several  examples  in  the  literature  of  single 
period  inventory  models  applied  to  military  logistics 
situations.  They  include  the  spare  parts  kit  problem  [Ref. 
7:pp.  281-295],  the  fly-away  kit  problem  [Ref.  6:p.  324], 
and  the  submarine  provisioning  problem  [Refs.  8:pp.  236-243, 
9:p.  407].  The  problem  with  each  of  these  military 
applications  is  that  they  assume  that  there  is  some  monetary 
stockout  cost  incurred  when  a  demand  exists  and  no  spares 
are  available  to  meet  the  demand.  This  is  realistic  in  the 
civilian  sector  where  real  monetary  values  can  often  be 
associated  with  each  stockout.  But,  in  most  practical 


military  applications  it  is  difficult,  if  not  impossible,  to 
determine  a  tangible,  realistic,  and  objective  stockout 
cost.  This  is  particularly  true  for  each  of  the  previously 
described  scenarios  because  there  is  no  way  to  place  a  cost 
on  the  inability  to  complete  a  mission  because  of  a 
deficiency  of  spares.  Thus,  it  is  difficult  to  apply  the 
models  existing  in  the  literature  to  these  real-world  Marine 
Corps  inventory  problems. 

This  problem  is  not  insurmountable.  In  the  civilian 
sector  when  an  organization  does  not  know  its  stockout  costs 
or  is  not  confident  in  its  estimates  of  them,  management 
will  set  service  levels  which  indicate  the  ability  to  meet 
customer  demands  from  stock  [Ref.  10:p.  149].  A  service 
level  is  basically  the  complement  of  the  probability  of 
stockout  and  really  represents  a  subjective  decision  by 
management  to  accept  some  degree  of  stockout  or  customer 
disservice  [Ref.  10:p.  150].  The  U.S.  Navy  actually  uses  a 
similar  concept  when  establishing  its  Aviation  Consolidated 
Allowance  Lists  ( AVCAL )  and  Coordinated  Shipboard  Allowance 
Lists  ( COSAL )  [Ref.  ll:pp.  4-30-4-50].  The  service  level 
measure  is  less  than  satisfactory  for  military  applications 
since  it  is  subjectively  established  by  combat  service 
support  planners.  Service  levels  do  not  provide  the 
operational  commander  in  the  field  with  a  meaningful 
indicator  of  the  impact  of  the  recommended  stockage  policy 


on  successful  mission  completion.  This  is  because  there  is 
no  direct  relationship  between  the  probability  of  stockout 
and  the  probability  of  successfully  completing  a  given 
mission.  A  high  service  level  does  not  equate  to  an  optimal 
probability  of  mission  completion.  This  is  particularly 
true  in  the  constrained,  multi-item  inventory  situation 
where  it  might  actually  be  better  to  accept  a  greater 
probability  of  stockout  ( lower  service  level  )  for  a 


secondary  reparable  with  a  low  failure  rate  in  order  to  be 
able  to  stock  more  of  another  secondary  reparable  with  a 
high  failure  rate. 

To  overcome  the  disadvantages  associated  with  the 
minimization  of  stockout  costs  (actually  variable  costs)  and 
service  level  approaches,  the  proposed  model  introduces  a 


new  measure  of  effectiveness,  the  probi 
completion  criven  the  specified  stockace 


of  mission 


i,  which  will 


be  referred  to  as  the  probability  of  mission  completion 
throughout  the  remainder  of  this  paper.  This  probability 
extends  the  familiar  concept  of  reliability  to  include  a 
system  and  its  inventory  of  spare  parts  [Ref.  8:p.  237].  A 
system  is  defined  to  be  the  quantity  of  a  principal  end  item 
required  to  be  operational  throughout  the  duration  of  a 


given  mission. 


This  end  item  must  be  a  "readiness 


reportable  pacing  item"  as  delineated  in  Marine  Corps 
Bulletin  3000  [Ref.  12:pp.  1-12  Enel  (2)].  A  readiness 


reportable  pacing  item  is  a  principal  end  item  whose 

operation  is  essential  to  the  ability  of  a  MAGTF  to  complete 

its  assigned  mission.  Simply  stated,  the  probability  of 

mission  completion  is  the  probability  that  k-out-of-n  (k<n) 

readiness  reportable  pacing  items  of  the  same  type  remain 

operational  for  a  specified  mission  with  a  given  stockage 

policy.  Black  and  Proschan  [Ref.  7:p.  283]  proposed  a 

similar  measure  of  effectiveness  in  their  early  work  using 

marginal  analysis  to  determine  optimal  spare  parts  kits  at 

minimum  cost.  They  called  it  probability  of  adequacy  or 

assurance  of  adequacy ,  which  is  defined  as 

The  probability  that  either  the  equipment  does  not  fail 
due  to  a.  .  .malfunction  or,  if  it  does  fail  for  this 
reason,  a  replacement  is  available  to  provide  continued 
successful  equipment  operation.  [Ref.  8:p.  237] 

However,  this  definition  applies  only  to  series  systems  in 

which  all  equipment  must  remain  operational  for  a  successful 

mission.  The  proposed  probability  of  mission  completion  is 

more  general  since  it  can  be  extended  to  "k-out-of-n 

systems".  This  is  required  in  the  present  study  because 

the  operational  commander,  not  the  military  logistician, 

determines  the  quantities  of  a  principal  end  item  which  he 

feels  are  necessary  to  complete  the  assigned  mission.  For 

example,  a  MAU  commander  may  feel  that  he  can  accomplish  a 

particular  mission  with  just  three  out  of  five  tanks  which 

permits  fewer  tank  secondary  reparables  to  be  stocked  than 

if  he  required  all  five  tanks  to  complete  the  same  mission. 


In  this  study  the  inventory  problems  in  each  of  the 
previously  described  scenarios  are  modeled  as  single  period, 
constrained,  multi-item  inventory  problems  in  which  the 
demands  are  the  result  of  independent  failures  of  secondary 
reparables.  In  the  mathematical  formulation  of  these 
models,  the  probability  of  mission  completion  is  maximized, 
as  opposed  to  the  standard  approach  of  minimizing  expected 
(variable)  costs.  This  provides  models  which  reflect  the 
fact  that  in  most  military  situations,  and  specifically  in 
the  three  Marine  Corps  scenarios  studied  here,  stockout  or 
shortage  costs  cannot  be  determined.  The  inclusion  and 
numbers  of  constraints  are  dependent  on  the  mode  of 
deployment.  However,  in  all  cases  studied,  they  are  linear 
and  will  be  referred  to  as  capacity  constraints.  Thus,  each 
model  will  be  viewed  as  a  mathematical  program  in  which  the 
objective  function  to  be  maximized  is  a  nonlinear, 
nondecreasing  function  of  the  decision  variable  representing 
the  quantity  of  secondary  reparables  to  stock.  Inclusion  of 
the  linear  capacity  constraint  or  constraints  result  in  an 
overall  integer  nonlinear  programming  problem.  Since  the 
objective  function  is  itself  a  complicated  function  of  the 
decision  variable,  the  solution  approach  is  to  linearize  the 
objective  function  and  then  reformulate  the  entire  problem 
using  binary  variables.  This  reformulation  converts  the 
original  integer  nonlinear  programming  problem  into  a  binary 


linear  programming  problem  which  is  then  coded  for  computer 
processing  using  the  Generalized  Algebraic  Modeling  System 
(GAMS)  LRefs.  13,  14,  15]  and  solved  using  the  Zero/One 
Optimization  Methods  (ZOOM)  [Ref.  16]  mixed  integer  program 
solver. 


Since  one  of  the  primary  purposes  of  this  study  is  to 
provide  the  combat  service  support  planners  of  the  Second 
Force  Service  Support  Group  at  Camp  Lejeune,  North  Carolina 
with  models  solvable  on  a  personal  computer,  it  is  desired 
to  avoid  dependency  on  commercial  optimization  packages  not 
locally  available.  Therefore,  the  initial  solution 
methodology  utilized  Langrangian  relaxation  to  establish  an 
upper  bound  on  the  optimal  solution  of  the  integer  linear 
problem.  Then  a  heuristic  based  on  marginal  analysis  is 
used  to  improve  upon  the  best  feasible  solution  obtained 


associated  with  the  search  procedure  used  in  establishing 
the  best  upper  bound  increases,  and  for  any  more  than  two 
constraints  actually  becomes  prohibitive.  So,  to  keep  the 
solution  approach  as  flexible  and  general  as  possible,  it 
was  decided  to  solve  all  integer  linear  problems  with  the 
GAMS/ZOOM  optimization  package  ( commerically  available  in 
March  1988),  thus  allowing  the  solution  of  models  with  more 
than  one  constraint  and  avoiding  restrictive  programming 
considerations . 

The  above  solution  approach  overcomes  the  disadvantages 
associated  with  more  standard  solutions  of  the  single  period 
inventory  model.  For  instance,  it  could  be  argued  that 
since  these  inventory  problems  are  all  examples  of  decision 
making  under  risk,  they  could  be  solved  using  decision 
theory  and  payoff  matrices  to  minimize  expected  costs. 
However,  as  was  pointed  out  earlier,  realistic  monetary 
stockout  costs  are  extremely  difficult  to  determine  in  most 
military  applications.  Additionally,  in  the  multiple  item 
inventory  problem,  the  payoff  matrices  become  very  large  and 
the  interrelationship  among  the  payoffs  associated  with 
different  combinations  of  items  becomes  difficult  to 
evaluate  and  enumerate,  so  it  is  computationally  more 
efficient  to  use  a  marginal  analysis  approach  [Ref.  10:p. 
269]  .  Marginal  analysis  provides  the  basis  for  the 
techniques  used  in  the  spare  parts  kit  problem  [Ref. 

14 


7 : pp . 283-284]  and  in  the  submarine  provisioning  problem 
[Ref.  8:p.  237],  but  it  has  practical  limitations  when  more 
than  one  constraint  must  be  considered.  Hadley  and  Whitin 
[Ref.  6 :  pp  328-330]  obtained  exact  integer  solutions  to 
their  single  constraint,  fly-away  kit  problem  using  dynamic 
programming  and  Dreyfus  and  Law  [Ref.  18:pp.  107-118]  offer 
dynamic  programming  algorithms  to  solve  single  constraint, 
cargo-loading  problems.  However,  they  all  point  out  that 
for  problems  involving  more  than  two  constraints  even 
dynamic  programming  has  computational  limitations, 
particularly  when  both  the  number  of  constraints  and 
decision  variables  are  large  [Refs.  6:p.  331,  18:pp*  107- 
118].  The  solution  approach  utilized  in  this  study  avoids 
these  shortcomings. 


C.  THESIS  OUTLINE 

This  thesis  develops  and  presents  a  method  for 
determining  spare  stockage  levels  of  secondary  reparables 
associated  with  the  readiness  reportable  principal  end  items 
of  a  given  Marine  Air-Ground  Task  Force.  In  Chapter  II,  the 
maintenance  data  base,  data  analysis,  and  problems 
encountered  during  this  analysis  are  discussed.  The 
objective  functions  are  defined,  developed,  and  justified  in 
Chapter  III.  Additionally,  each  scenario  is  modeled  and 
formulated  as  an  integer  nonlinear  programming  problem  ror 


which  the  objective  function  is  linearized  so  that  the 
problem  can  be  reformulated  as  a  integer  linear  problem  to 
facilitate  its  solution.  Specific  solution  methodology  is 
presented  in  Chapter  IV.  The  initial  approach  using 
Lagrangian  relaxation  with  a  heuristic  is  briefly  described 
before  the  more  general  procedure  using  GAMS/ZOOM  is 
developed  and  implemented.  In  Chapter  V  results  using 
actual  data  for  a  typical  Marine  Amphibious  Unit  are 
presented  and  verification  of  the  model  is  conducted  using 
different  sets  of  simulated  failure  data.  The  chapter 
concludes  with  a  discussion  of  computational  experience  for 
both  the  mainframe  and  personal  computer.  Chapter  VI 
contains  conclusions  and  recommendations.  Finally,  listings 
of  both  FORTRAN  77  and  GAMS  source  codes  are  provided  in  the 
appendices . 


II.  DATA  ISSUES  AND  ANALYSIS 


A  primary  impetus  for  this  study  was  the  desire  to 
develop  a  reliability-based,  spares  stockage  model  for 
secondary  reparables  utilizing  the  maintenance  data  resident 
in  the  Marine  Corps  Integrated  Maintenance  Management  System 
( MIMMS ) .  Currently,  spare  stockage  models  within  the  Marine 
Corps  are  limited  to  those  based  on  forecasts  from  past 
demands  generated  against  the  supply  system.  This 
forecasting  approach  neglects  to  consider  that  the  source  of 
the  demands  is  the  failure  of  equipment  such  as  secondary 
reparables,  and  it  is  precisely  this  type  of  failure 
information  that  can  be  extracted  from  relevant  MIMMS  sub¬ 
files  which  permits  the  development  of  a  reliability-based 
model . 

A.  DESCRIPTION  OF  DATA 

MIMMS  is  an  automated  system  which  accumulates 
maintenance  data  on  all  serialized,  principal  end  items  in 
the  Marine  Corps  inventory.  It  records  and  accumulates,  in 
various  sub-files,  maintenance  actions  performed  on  a 
serialized  end  item  during  its  service  life  at  each  of  the 
three  levels  of  maintenance  (i.e.,  organizational, 
intermediate,  and  depot)  [Ref.  19].  In  particular,  the 
MIMMS  sub-file  of  relevance  to  this  study  is  the  Equipment 


Repair  Order  History  File.  This  sub-file  maintains  all 
equipment  repair  orders  ( EROs )  that  were  completed  on  every 
serialized  end  item  during  the  past  eighteen  months,  and  it 
is  updated  on  a  quarterly  basis.  An  ERO  is  opened  each  time 
a  serialized  end  item  is  inducted  into  the  maintenance  cycle 
and  contains  information  such  as  status  of  the  end  item 
(i.e.,  combat  deadlined,  safety  dendlined,  mission  degraded, 
etc.  ),  the  meter  reading  at  the  time  of  the  current 
maintenance  action  (expressed  in  either  hours,  miles,  rounds 
or  days  depending  on  the  end  item),  and  the  defect  code 
which  describes  the  nature  of  the  failure.  By  looking  at 
only  those  EROs  in  the  MIMMS  ERO  History  File  which  were 
opened  for  critical,  combat  deadlining  repairs  (i.e., 
category  code  M),  a  new,  local  user  file  can  be  created. 
This  user  file  is  organized  by  principal  end  item  and 
includes  only  those  end  items  which  are  considered 
"readiness  reportable"  [Ref.  12:pp.  Enel  (1)  1-23). 
Associated  with  each  end  item  are  all  serial  numbers  for 
which  a  combat  deadlining  ERO  exists;  each  serial  number  has 
a  meter  reading  and  defect  code  to  reflect  each  ERO 
completed  on  that  particular  serial  number  during  the 
eighteen  month  period  covered  by  the  ERO  History  File. 
Obviously,  there  may  be  multiple  entries  under  the  same 
serial  number  indicating  that  this  end  item  suffered  more 
than  one  combat  deadlining  failure  during  the  reporting 


experienced  by  the  particular  serial  number.  Table  1 
provides  a  sample  extract  from  this  user-created  data  file 
and  displays  its  contents  for  several  tank  serial  numbers. 


TABLE  1 

EXTRACT  FROM  LOCAL  USER  DATA  FILE 


SERIAL  # 

EOTC* 

METER 

READING 

DEFECT 

CODE** 

502814 

H 

74 

D05 

(POWER  TRAIN, CLUTCH) 

502828 

H 

51 

A34 

(ENGINE  REPLACE) 

H 

75 

J34 

(COOLING  SYS  REPLACE) 

H 

107 

C34 

(POWER  RACK  REPLACE) 

H 

155 

A34 

(ENGINE  REPLACE) 

H 

163 

B34 

(TRANSMISSION  REPLACE) 

502831 

H 

81 

A34 

(ENGINE  REPLACE) 

H 

149 

B34 

(TRANSMISSION  REPLACE) 

502835 

H 

62 

KOI 

(ELEC  SYS,  GENERATOR) 

*  EOTC=  Equipment  Operating  Time  Code 

(H=Hours,  M=Miles,  R=Rounds,  D=Days) 

**  Defect  Codes  are  alphanumeric  codes  and  have  been 
annotated  here  for  display  purposes  only. 


The  data  in  this  consolidated  file  can  be  used  to  estimate 
failure  rates  for  each  secondary  reparable  in  either  of  two 
ways  and  both  methods  will  be  described  in  the  f  -lowing 
section  along  with  their  advantages  and  disadvantages. 
It  must  be  stated  that  both  the  format  and  method  of  update 
for  many  of  the  MIMMS  files  require  improvement  before  a 
truly  comprehensive  analysis  of  the  data  is  possible.  A 
case  in  point  is  the  ERO  History  File  used  in  this  study, 
and  recommendations  for  its  improvement  are  provided  in  the 


following  paragraph  and  in  Chapter  VI  which  minimize  the 
impact  on  current  MIMMS  file  management  and  the  amount  of 
effort  required  at  the  using  unit  input  level. 

One  of  the  basic  problems  with  the  ERO  History  File  is 
that  potentially  useful  data  is  being  discarded  during  each 
quarterly  update  in  order  to  minimize  file  size.  What  is 

proposed  here  is  to  maintain  a  local  (i.e.,  at  the  Force 

Service  Support  Group  level)  consolidated  version  of  the  ERO 
History  File  with  a  format  similar  to  that  described  in  the 
preceding  paragraph.  This  file  should  be  updated  when  its 
parent  ERO  History  File  is  updated  at  the  end  of  each 
quarter.  Since  this  consolidated  file  contains  only 

information  extracted  from  combat  deadlining  EROs,  which 
comprise  a  small  fraction  of  the  total  EROs  in  the  parent 
file,  there  is  a  greater  capability  to  accumulate  relevant 
data  without  the  need  to  automatically  discard  the  oldest 
three  months  of  data.  The  basic  algorithm  for  each 

quarterly  update  is: 

Step  1.  Select  a  combat  deadlining  ERO  (Category  Code  M) 
from  the  ERO  History  File  if  one  exists.  If  none 
exist,  STOP. 

Step  2.  Determine  if  the  serial  number  associated  with 
this  ERO  is  already  resident  on  the  consolidated 
file.  If  it  is,  continue  with  Step  3.  If  it  is 
not,  add  the  serial  number  and  its  associated 
meter  reading  and  defect  code  to  the  consolidated 
file  under  the  appropriate  principal  end  item. 
Return  to  Step  1. 


Step  3.  Determine  if  the  same  defect  code  exists  for  this 
serial  number.  If  it  does,  continue  with  Step  4. 
If  it  does  not,  add  the  meter  reading  and  defect 
code  to  the  list  associated  with  this  serial 
number.  Return  to  Step  1. 

Step  4.  Subtract  the  two  meter  readings  to  obtain  a  time 
between  failure  for  this  type  of  defect.  Add 
this  new  data  point  to  the  list  associated  with 
this  type  of  defect.  Delete  the  record 
associated  with  the  older  meter  reading  (i.e., 
the  smallest  meter  reading)  and  add  the  new  meter 
reading  and  defect  code  to  the  end  of  the  list 
associated  with  this  serial  number.  Return  to 
Step  1. 

The  result  is  an  updated  version  of  the  consolidated  file 
which  is  retained  at  the  local  level  for  update  during  the 


next  quarterly  cycle.  Each  readiness  reportable  principal 
end  item  in  the  file  is  broken  down  by  serial  number  which 
in  turn  contain  records  of  the  most  recent  meter  readings 
and  defect  codes.  Each  end  item  also  consists  of  lists 
containing  times  between  failures  for  each  type  of  defect. 
Thus  this  local  consolidated  file  does  not  alter  any 
existing  MIMMS  file  definitions,  nor  does  it  necessitate  the 
creation  of  any  new  system  file  definitions.  More 
importantly,  it  provides  the  combat  service  support  planner 
with  a  consolidation  of  all  the  pertinent  information 
required  to  compute  failure  rates  and  track  the  life  history 
of  any  secondary  reparable  associated  with  a  readiness 
reportable  end  item.  This  is  a  significant  improvement  over 
the  present  method  which  is  restricted  to  using  the  most 
recent  eighteen  months  of  failure  data  resident  in  the  ERO 


History  File,  and  it  is  this  restriction  that  limits  the 
accumulation  of  failure  times  for  certain  types  of  secondary 
reparables  (e.g.,  certain  radio  transmitters). 

B.  DATA  ANALYSIS 

Using  the  data  resident  in  the  consolidated  file,  the 
operating  times  between  failures  for  each  type  of  secondary 
reparable  were  computed.  Failure  rates  were  estimated  using 
two  different  approaches.  In  the  first  approach,  an 
exponential  lifetime  distribution  was  assumed  for  each 
secondary  reparable  and  the  maximum  likelihood  estimator 
(M.L.E. )  for  the  parameter  A,  of  an  exponential  distribution 
was  used  to  provide  an  estimate  of  the  failure  rate.  The 
second  approach  involves  a  Bayesian  analysis  of  the  times 
between  failure  to  determine  a  posterior  distribution  for 
the  failure  rate  of  a  secondary  reparable. 

In  deriving  the  M.L.E.  for  A^,  it  is  assumed  that  the 
observed  lifetimes  of  each  secondary  reparable  i  form  a 
relevant  random  sample  of  size  nj_  from  an  exponential 
distribution  with  unknown  parameter  /^ .  The  likelihood 
function  must  also  account  for  the  exposure  times  during 
which  no  failures  of  type  i  occurred.  If  these  terms  are 
not  considered,  the  computed  failure  rate  will  tend  to  be 
too  large  and  excessive  stockage  levels  would  result.  For 
the  n.^  observed  failure  times  and  the  observed  exposure  time 


t'i  during  which  no  failures  occurred,  the  likelihood 
function  is: 
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Maximizing  the  likelihood  function  with  respect 

yields  the  following  estimator  of  the  failure  rate: 

a  n,  rtj 
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where  represents  the  total  exposure  time  for  secondary 
reparable  i.  [Ref.  20 : pp. 282-296] 

The  Bayesian  approach  also  assumes  that  the  distribution 
of  the  lifetimes  for  each  secondary  reparable  i  is 
exponential  with  parameter  /  i;  however,  it  further  assumes 
that  the  exact  value  of  this  parameter  is  a  realization  of  a 
random  variable.  Even  though  ).  ^  is  not  precisely  known,  it 
has  a  prior  distribution  which  is  taken  here  to  be  a 
conjugate  Gamma  distribution  with  shape  parameter 
a  i(ai>0)  and  scale  parameter  /?i(  />i> 0): 
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The  prior  distribution  parameters  are  estimated  using  a 
parametric  empirical  Bayes  ( PEB )  approach  for  failure  rates 
as  outlined  by  Gaver  and  Lehoczky  [Ref.21:pp.  220-224].  The 


PEB  approach  uses  the  entire  data  set  to  compute  estimates 

A  A 

a  and  which  are  substituted  into  the  applicable 

formulas.  To  compute  the  prior  parameters  in  PEB  requires 
the  maximization  of  the  "marginal  likelihood  function": 
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where, 

m  individual  source  of  data  for  component  i  (e.g., 
tank  serial  #502826) 

A/,  total  number  of  data  sources  for  component  i 

/-im  failure  rate  for  component  i  of  data  source  m 

hm  observed  exposure  time  of  component  i  on  data  source 
m 

nim  number  of  failures  of  component  i  on  data  source  m 

ai  shape  parameter  of  Gamma  prior  distribution  of 

failure  rate  for  component  i 

Pi  scale  parameter  of  Gamma  prior  distribution  of 

failure  rate  for  component  i. 

Integration  results  in  the  following  marginal  likelihood 
function  for  secondary  reparable  i: 


A U 

*  Pi  ■  %i  *  t-i )  =  |  j 

m=  1 


+  rjP 

nim !  r  (*/) 


(2.4) 


which  is  the  product  of  individual 

Pi 

distributions  with  parameters  p  = - — 

lim  •  Pi 


negative 


binomial 


and 


r  =  «. 


The  numerical  maximization  of  this  function  subject  to  known 
upper  and  lower  bounds  on  each  a  ^  and  p  ^  is  discussed  in 
Chapter  IV.  The  important  result  of  this  maximization  is  an 

A  A 

estimate  for  the  shape  (  aj_)and  scale  (p±)  parameters 
of  the  Gamma  prior  distribution  that  incorporates  the  actual 
observed  operational  data  for  each  secondary  reparable. 

Returning  to  the  Bayesian  analysis,  the  ni  observed 
times  between  failure  and  the  observed  exposure  times 
during  which  no  failures  occurred  are  all  a  function  of 
/  and  give  rise  to  the  likelihood  function. 


f{u  I  '-*) 


.n: 


=  //exp 


(2.5) 


where  r(-  once  again  represents  the  total  exposure  time  for 
secondary  reparable  i.  We  are  now  interested  in  the 
updated  distribution  of  /.  i  after  the  times  between  failure 
have  been  observed  so  as  to  provide  a  current  estimate  of 
future  failure  events.  This  is  the  posterior  distribution 
and  it  is  proportional  to  the  product  of  the  likelihood 
function  (2.5)  and  the  prior  distribution  (2.2): 


fUi  I  Li)  /(£,  I  /)/(/)  ■ 


Thus, 


where  ci  is  a  constant  factor  which  depends  on  the  observed 
times  between  failure,  but  does  not  depend  on  A^.  The  value 
of  this  constant  can  be  determined  using  normalization.  The 
resulting  posterior  distribution  is: 


ft'-i  I  £,•)  =  /U/!  data)  = 


i  11 11 _  i  nt  +  m  —  1 

‘  ,  A  .  ' 


(t.  +  /?,)/, 


Tin,  4-  af) 


a  Gamma  distribution  with  shape  parameter  «,  +  a,-  and  scale 
parameter  r  ,•  +  /?,•  .  [Ref.  20 :  pp .  257-280] 

The  posterior  distribution  could  also  have  been  derived 
directly  using  the  fact  that  the  family  of  Gamma 
distributions  serves  as  a  conjugate  family  of  prior 
distributions  for  samples  from  an  exponential  distribution. 
Thus,  using  the  theorem  for  Gamma  conjugate  priors  [Ref. 
20:pp.  271-272],  the  posterior  distribution  for  / ^  given 

the  observed  times  between  failure,  is  a  Gamma  distribution 
with  shape  parameter  n^  +  and  scale  parameter  r(-  +  /?,•  . 

This  agrees  with  the  result  derived  in  the  previous 
paragraph  and  indicates  that  the  posterior  distribution  can 
be  defined  knowing  only  the  number  of  observations  and  the 
sum  of  all  the  exposure  times.  This  is  an  important  result 
because  it  enables  the  posterior  distribution  to  be  updated 
very  easily  as  new  observations  are  obtained. 

Since  the  posterior  distribution  of  A,-  has  been  shown 


to  be  a  Gamma  distribution,  its  expected  value  is: 


,  rl:  +  U: 

E [A,-  I  data]  =  — - -  . 

-  h 

This  expected  value  also  happens  to  be  a  "Bayes  estimator" 
of  /.  when  a  "squared  error  loss  function"  is  used  [Ref. 
20:p.  227],  This  estimator  could  also  have  been  used  to 

provide  an  estimate  for  each  /.  However,  the  primary 

purpose  of  the  Bayesian  analysis  was  not  to  provide  a  point 
estimate  of  but  rather  to  provide  a  distribution  (i.e., 

posterior  distribution)  for  each  /.  ^  which  could  be  used  to 
make  probabilistic  predictions  about  demand  functions 


III.  MODEL  FORMULATION  AND  DESCRIPTION 


In  this  chapter,  the  two  different  objective  functions 
considered  in  this  study  are  defined,  derived,  and 
justified.  Additionally,  each  of  the  three  scenarios 
described  in  the  introduction  is  formulated  as  an  integer 
nonlinear  program  whose  objective  function  is  subsequently 
linearized  to  create  an  integer  linear  problem. 

A.  MODEL  DEVELOPMENT 

The  following  formulation  is  developed  to  model  the 
determination  of  secondary  reparable  inventory  levels  for 
each  of  the  three  scenarios  described  in  the  introduction. 
With  the  exception  of  the  number  of  constraints,  the 
formulation  is  the  same  for  each  scenario. 

INDICES: 

/=  1 . A'  secondary  reparable  or  component 

DATA: 

W(  unit  shipping  weight  of  secondary  reparable  i 

[lbs] 

W  weight  capacity  of  aircraft  or  ship  [lbs] 

vi  unit  volume  of  secondary  reparable  [cuft] 

V  volume  capacity  of  aircraft  or  ship  [cuft] 

ui  maximum  allowable  number  of  spares  for 

component  i 


FUNCTIONS: 


pt{xj)  probability  that  the  number  of  readiness 

reportable  pacing  items  (specified  by  the 
operational  commander)  containing  secondary 
reparable  i  will  survive  a  mission  given  xi 
spares  are  stocked 


DECISION  VARIABLE: 

Xj  Number  of  spares  to  stock  for  secondary 

reparable  i 


FORMULATION: 

v 

max  (pi) 

i- 1 
.V 

s.t.  yw^iv 

i=i 

v 

/=  1 

o  ^  <  “i  /  =  1, ...  ,.v 

xi  integer 


The  objective  function  is  a  mathematical  expression  of 
the  probability  that  a  quantity  (specified  by  the 
operational  commander)  of  each  readiness  reportable  pacing 
item  in  a  MAGTF's  Table  of  Equipment  will  survive  a  mission 
of  specified  duration  given  a  particular  stockage  policy  for 
secondary  reparables.  This  probability  was  referred  to  in 
the  introduction  as  the  probability  of  mission  completion. 


<  2 

£;  l 

rj  z 


In  terms  of  this  definition,  successful  mission  completion 
hinges  on  the  survival  of  the  MAGTF's  pacing  items.  Each 
pacing  item  is  composed  of  secondary  reparables.  The 
failure  of  any  one  secondary  reparable  will  cause  the  pacing 
item  to  be  combat  deadlined  and  unable  to  complete  its 
mission.  Thus,  each  pacing  item  is  viewed  as  a  system  of 
secondary  reparables  in  series.  However,  in  a  MAGTF's  Table 
of  Equipment,  there  is  not  just  one  of  each  type  of  pacing 
item,  but  rather  a  specific  quantity  of  each  type  necessary 
to  accomplish  the  assigned  mission.  This  minimum  quantity 
required  will  not  only  vary  from  one  mission  to  another,  but 
also  among  different  MAGTF  commanders  when  confronted  with 
identical  missions.  Therefore,  the  different  pacing  items 
taken  as  a  whole  do  not  constitute  a  simple  series  system  in 
which  all  components  must  survive  the  given  mission. 

In  defining  the  objective  function,  two  different 
functional  forms  were  derived  and  studied.  Both  rely  on  the 
assumptions  that  secondary  reparables  fail  at  constant  rates 
and  that  the  failures  of  individual  secondary  reparables 
are  independent  of  one  another.  The  independence  assumption 
is  certainly  true  for  the  secondary  reparables  associated 
with  different  types  of  pacing  items;  however,  it  may  not 
necessarily  be  true  for  the  secondary  reparables  of  a 
specific  pacing  item.  For  example,  the  failure  of  M60A1 
tank  engines  is  independent  of  the  failure  of  LVTP7A1 


amtrack  engines,  but  the  failure  of  amtrack  final  drives  may 
not  be  independent  of  amtrack  transmission  failures.  These 
assumptions  together  with  the  fact  that  secondary  reparables 
from  one  type  of  pacing  item  are  not  interchangeable  with 
those  of  another  (e.g.,  a  tank  engine  cannot  be  substituted 
for  an  amtrack  engine)  make  it  possible  to  consider  only 
secondary  reparables  without  having  to  specify  the  pacing 
items  of  which  they  are  components. 

The  functional  form  of  the  first  objective  function,  in 
which  the  failure  rate  of  each  secondary  reparable  is 
assumed  to  be  known,  is  as  follows: 


ai  +  */  “  ki 


>»(*»)  =  X 


(aj.jtjf 


cxp-^'' 


where , 


number  of  pacing  items  containing  secondary  reparable 


ki  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (ki  <  a,) 

A 

/,■  M.L.E.  of  failure  rate  for  secondary  reparable  i 

f(  mission  duration  during  which  secondary  reparable  i 
must  operate. 

This  function  is  a  direct  consequence  of  the  constant 
failure  rate  assumption  for  each  secondary  reparable. 
Underlying  this  constant  failure  rate  is  an  inherent 


exponential  lifetime  distribution  which  implies  that  the 


number  of  failures  defines  a  Poisson  (arrival)  process  with 


exponential  interarrival  times.  The  number  of  failures  of 
secondary  reparable  i  to  occur  in  any  fixed  interval  of 
time,  such  as  a  mission  of  duration  t±,  will  have  a  Poisson 
distribution  with  mean  /  However,  since  there  are  a 
total  of  a_£  pacing  items  containing  secondary  reparable  i 
and  operating  during  this  mission,  the  total  number  of 
failures  will  be  the  sum  of  independent,  Poisson  random 
variables  (each  with  a  mean  ±tj_)  which  is  itself  a  Poisson 
random  variable  with  a  mean  .  This  is  really  only 
an  approximation  since  not  all  a^  pacing  items  may  be 
required  during  the  mission.  In  which  case  the  overall 
failure  rate  would  be  smaller  implying  that  we  are  probably 
overestimating  the  number  of  spares  required.  On  the  other 
hand,  it  can  be  argued  that  with  fewer  than  a items  the 
overall  failure  rate  of  the  operational  items  may  increase 
since  the  same  mission  must  be  accomplished  with  fewer 
assets.  In  any  case,  this  can  be  viewed  as  a  mild,  worst 
case  approximation. 

Therefore,  equation  (3.1)  represents  the  probability 
that  there  are  less  than  or  equal  to  some  specified  number 
of  failures  of  secondary  reparable  i  during  a  mission  of 
duration  t^.  This  specified  number  is  the  sum  of  the  number 
of  spares,  and  the  difference  between  the  number  of 
pacing  items  containing  secondary  reparable  i  (a^)  and  the 


required  minimum  number  of  pacing  items  containing  secondary 
reparable  i  ( k^_ ) .  As  an  example,  consider  the  situation 
where  the  Table  of  Equipment  of  some  MAGTF  consists  of  five 
tanks  each  containing  one  engine  (  a^  = 5)  and  assume  that  the 
MAGTF  commander  decides  that  to  accomplish  the  assigned 
mission  he  must  have  a  minimum  of  four  tanks  operational  at 
all  times  (k^= 4).  Then,  assuming  one  spare  tank  engine 
(Xj_  =  l)  is  stocked,  there  can  be  no  more  than  two  engine 
failures  ( a x ^-k 5+l-4=2)  to  adhere  to  the  commander's 
guidance.  It  can  be  shown  that  this  probability  is  equal  to 
the  probability  that  the  quantity  of  each  type  of  pacing 
item  (or  equivalently  all  its  component  secondary 
reparables)  remains  above  a  specified  level  (ki),  because 
the  following  relation  must  always  hold  true: 

(number  of  failures  of  secondary  reparable  i}  + 

(number  of  secondary  reparables  i  operational  and  in 
inventory}  = 

(total  number  of  secondary  reparables  i,  ( a^  +  x^ ) } 

Thus, 

P( (number  of  secondary  reparables  i  operational  and 
in  inventory  at  t }>k.j_)  = 

?((a^+Xy}  -  (number  of  failures  of  secondary  reparable 
i  in  z)>k;i_)  = 

P( (number  of  failures  of  secondary  reparable  i  in  t} 

iai+xi“ki ) • 

Therefore,  the  probability  of  mission  completion  for  the 
entire  MAGTF  is  the  product  of  the  probabilities  that  the 


individual  secondary  reparables  exceed  a  level,  k^,  given  x^ 
spares  (or  equivalently  experience  at  most  a^+x^-k^ 
failures).  So,  the  maximum  number  of  operational  pacing 
items  of  a  specific  type  at  any  given  moment  will  be  equal 
to  the  lowest  level  of  its  component  secondary  reparables. 
It  is  this  minimum  level  that  must  exceed  the  specified 
level  k since  if  the  minimum  level  exceeds  k^r  all  levels 
will  exceed  Now,  by  the  assumption  that  the  secondary 

reparable  failures  are  independent  combined  with  the  fact 
that  if  one  type  of  pacing  item  cannot  complete  the  mission 
the  entire  MAGTF  will  not  be  able  to  complete  the  mission, 
the  probability  that  the  minimum  level  exceeds  the  specified 
level  (i.e.,  Pj_(x±))  can  be  expressed  as, 

.v 

Pi  (.tj )  •  p2(x2) .....  Px(x y)  =  j  3  2) 

1=1 

which  accounts  for  the  functional  form  of  the  objective 
function. 

Similar  to  the  first  objective  function,  the  second 
objective  function  also  takes  advantage  of  the  assumptions 
of  constant  failure  rates  for  each  secondary  reparable  and 
independence  among  the  failures  of  different  secondary 
reparables;  however,  it  incorporates  a  more  realistic 
approach  as  to  the  nature  of  the  failure  rates.  In  the 
previous  objective  function,  the  failure  rate  A^,  of  each 


secondary  reparable  was  assumed  to  be  known  (actually  an 


estimate  computed  using  the  maximum  likelihood  estimator  of 
the  failure  rate  must  realistically  be  used).  For  this 
second  objective  function,  a  Bayesian  approach  is  taken  in 
which  the  exact  value  of  ^  is  explicitly  recognized  to  be 
unknown.  The  observed  data  (i.e.,  times  to  failure)  are 
used  to  update  the  prior  information  concerning  each  failure 
rate  to  form  a  posterior  distribution  as  discussed  in  the 
chapter  on  data  analysis.  This  posterior  distribution  is 
necessary  to  derive  the  predictive  distribution  and  the 
functional  form  of  p^(x^),  since  we  are  interested  in 
computing  the  probabilty  distribution  of  the  number  of 
failures  of  secondary  reparable  i  which  occur  during  some 


mission  of  length  t^.  By  analogy  to  the  previous  derivation 
of  ,  it  is  desired  to  consider  this  number  of 

failures  because  we  want  this  quantity  to  be  less  than  some 
specified  level  (i.e.,  a^+x^-k^)  so  as  to  insure  the  entire 
MAGTF  completes  the  assigned  mission.  Assuming  that  the 
failure  rate  /.  ^  is  given,  the  total  number  of  failures  of 
aj_  secondary  reparables  of  type  i  in  a  mission  of  length  tj_ 
will  have  a  Poisson  distribution  with  approximate  mean 

ai  '■  itir 

/’(Exactly  fx  failures  in  = 

But  /  ^  is  not  precisely  known  and  has  the  following 
posterior  density  (as  derived  in  the  data  analysis  chapter) 
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given  the  observed  times  to  failure  and  prior  information: 


C  -»■ .  _i_  /?  )ni  a  A 

JUi  I  data)  =  — - — - +  *«'  - ' 1  cxp_/-'(T'  +  ^ 

r>,  + «,) 


where. 


v  total  exposure  time  for  secondary  reparable  i 


ni  number  of  observations  for  secondary  reparable 


a,-  empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(  A,  :  a,  ,  Pi  )  prior  distribution 


;3{-  empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(  A, ;  a, ,  )  prior  distribution  . 


This  can  be  rewritten  as. 


/(A,- 1  data)  = 


+  ft))”' +  ‘‘  _1(T,  +  1)  exp-;-><Ti +  « 


r(«i  +  a,) 


Now,  to  calculate  the  predictive  distribution, 


/’(Exactly^  failures  in  (0,rj)  =  f  /’(Exactly  fi  failures  in  (O.rJ  I  A,)  J[/.t  I  data WA, 

*>o 


which  is  upon  substitution. 


LL 

.n  '  n 


.  -a/VA/  111  t-fftLl - 1  t  exn-^<Ti+  J/) 


n«/  +  a,) 


Evaluation  of  the  integral  yields: 


/■(Exaalvy;  failures  in  (0,,,])  =  +  *'  ( - 2£ - f 

V<V|.+  ri  +  Jiy  a/,  +  T,  +  ?,  ’ 


This  is  a  negative  binomial  distribution  with  parameters 


ri +  Pi 


aih  "E  T<  "E  Pi 


and  r  =  nt  +  a . 


r«.»y,o 


a 
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Even  though  p^(x^)  is  defined  as  the  probability  that 
the  stockage  of  x^  spares  for  secondary  reparable  i  allows 
the  pacing  item  of  which  it  is  a  component  to  exceed  some 
minimum  level  ( k± ) ,  it  has  been  shown  that  this  is 
equivalent  to  the  probability  that  the  number  of  failures  of 
secondary  reparable  i  during  mission  t ^  is  less  than  or 
equal  to  a^+x^-k^  (i.e.,  P( number  of  failures  of  secondary 
reparable  i  in  (0,tj_]<  a^+x^-k^)).  Therefore,  the  second 
functional  form  of  p^(x^)  is: 

aiti 


<*ih  +  Ti  +  Pi 

where, 

ai  number  of  pacing  items  containing  secondary  reparable 
i 

ki  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (k^<a^) 

number  of  observations  ( data  points )  of  times  between 
failure  for  secondary  reparable  i 

Ti  total  observed  exposure  time  for  secondary  reparable  i 

^  mission  duration  during  which  secondary  reparable  i 
must  operate 

at-  empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(  A, ;  a,- ,  Pi  )  prior  distribution 

A 

Pi  empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(  A( ;  a,- ,  Pi  )  prior  distribution  . 


a -t-  x:  —  ki 


fa)  =  y 


r  (rt;  +  y-i  +  tn) 


+  Pi 


w=o  tnl  r(«(-  -  a() 


aih  +  T«  +  Pi 


r'( 


The  reasoning  behind  the  functional  form  of  equation  (3.2) 
developed  for  the  first  pi(xi)  still  holds  true  for  the 


newly  derived  p^(x^);  therefore,  the  functional  form  of  this 
objective  function  is  still  the  product  of  the  individual 
' s. 

The  linear  weight  capacity  restricts  the  total  shipping 
weight  (i.e.,  secondary  reparable  plus  container)  for  all 
secondary  reparables  stocked  to  be  less  than  or  equal  to  the 
weight  capacity  of  the  ship  or  airplane.  Similarly , the 
linear  constraint  on  volume  or  space  capacity  requires  that 
all  secondary  reparables  stocked  not  exceed  the  volume 
available.  The  final  constraint  in  PI  indicates  that  there 
is  an  upper  bound  on  the  number  of  spares  for  each  secondary 
reparable  that  can  be  stocked.  This  is  a  very  real 
constraint  because  there  are  only  a  finite  number  of  each 
type  of  secondary  reparable  at  the  wholesale  inventory  level 
from  which  the  retail  inventories  used  to  support  all 
exercises,  deployments,  and  contingencies  can  be 
established.  Thus,  a  combat  service  support  planner  must 
consider  the  requirements  of  the  force  as  a  whole  so  as  not 
to  commit  too  many  secondary  reparables  to  any  one 
exercise . 

The  number  of  constraints  is  dependent  on  the  scenario. 
Any  scenario  that  involves  the  airlift  of  the  ground  combat 
element  and  its  supporting  forces  into  the  exercise  or 
operation  area,  such  as  the  short-term  exercises  described 
in  the  introduction,  will  have  only  one  constraint.  This  is 


the  constraint  on  the  allowable  cargo  load  (ACL).  Allowable 
cargo  load  is: 

The  initial  limiting  factor  in  determining  airlift 
requirements.  .  .  .  This  is  a  weight  limiting  factor 

expressed  in  terms  of  STON  [short  tons]  which  cannot  be 
exceeded  in  planning  aircraft  loads.  The  ACL  is 
established  by  the  Air  Force  for  each  type  aircraft  for  a 
specific  mission.  [Ref.  22:pp.  4-9] 

Since  the  formulation  of  problem  PI  now  reduces  to  just  one 

constraint,  it  is  recognizable  as  a  nonlinear  knapsack 

problem  due  to  the  nonlinear  objective  function.  The 

individual  secondary  reparables  represent  the  items  or 

commodities  that  must  fit  into  a  knapsack  while  the  ACL 

constraint  corresponds  to  the  knapsack's  weight  or  volume 

capacity. 

The  other  two  scenarios  described  in  the  introduction 
both  involve  two  constraints.  These  two  constraints 
represent  the  weight  and  volume  capacities  of  the  shipping 
assigned  to  perform  the  sealift.  In  the  case  of  the  MAUs 
deploying  to  the  Mediterranean  Sea,  the  secondary  reparables 
are  normally  stored  aboard  one  of  the  ships  of  the 
amphibious  task  force  such  as  a  LHA,  LPH,  LKA,  LSD,  or  LST. 
Thus,  the  constraints  reflect  the  weight  and  volume  capacity 
of  the  ship's  cargo  hold  or  compartment  which  has  been 
designated  for  secondary  reparable  storage.  On  the  other 
hand,  in  considering  secondary  reparable  support  of  the  MPF 
MAB,  the  secondary  reparables  are  spread- loaded  among 
several  of  the  MPF  squadron  ships;  therefore,  the  weight  and 


volume  constraints  are  actually  the  combined  weight  and 
volume  capacities  of  the  cargo  holds  of  the  specified  MPF 
squadron  ships.  Since  both  of  these  sealift  scenarios 
involve  two  contraints,  problem  PI  can  no  longer  be  regarded 
as  a  nonlinear  knapsack  problem,  but  it  is  still  an  integer 
nonlinear  programming  problem. 

Regardless  of  the  number  of  constraints  or  whether 
equation  (3.1)  or  (3.3)  is  used  as  the  functional  form  for 
p±(Xj_),  the  objective  function  is  a  complex,  nonlinear 
function  of  the  decision  variable.  For  this  reason,  problem 
PI  will  be  transformed  into  a  binary  linear  programming 
problem.  The  following  section  explains  this  transformation 
and  the  new  formulation. 

B.  TRANSFORMATION  TO  A  BINARY  LINEAR  PROGRAM 

It  can  be  seen  from  the  final  constraint  of  problem  PI 
that  the  decision  variable,  x^,  is  constrained  to  a  set  of 
discrete  values.  Therefore  the  objective  function  will 
first  be  separated  using  the  properties  of  logarithms  and 
then  the  transformation  technique  for  converting  discrete 
variables  to  binary  variables  will  be  used  to  reformulate 
the  problem  [Ref.  23:p.  12].  Taking  the  natural  logarithm 
of  the  objective  function  yields: 


Now,  the  final  constraint  in  PI  can  be  expressed  as. 


Xj  £  S;  {  SjQ  >  -5;  1 1 


•S/u..}  —  {0,1,  ,Wj} 


and  this  is  equivalent  to  the  constraint  set  [Ref.  23  :p. 
12]  : 
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Applying  the  separation  technique  and  transformation  to 
problem  PI  results  in  the  new  formulation: 


INDICES: 

/=  1....  ,.V 

7  =  0 . «i 

DATA: 

h 


secondary  reparable  or  component 
number  of  secondary  reparables  i 

ln( P±( j  ) ) 

total  shipping  weight  for  j  spares  of  type  i 
[lbs] 


W 

vij 

V 

ui 


weight  capacity  of  airplane  or  ship  [lbs] 

total  volume  of  j  spare  of  typei[cuft] 

volume  capacity  of  airplane  or  ship  [cuft] 

maximum  allowable  number  of  spares  for 
component  i 


DECISION  VARIABLE: 


*  -  {J 


if  j  spares  of  type  i  are  selected 
otherwise 


FORMULATION: 
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Assuming  known  failure  rates,  the  coefficients,  in  the 

objective  function  are: 
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where, 

ai  number  of  pacing  items  containing  secondary  reparable 
i 

kj  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  ( k^£a^ ) 

A 

M.L.E.  of  failure  rate  for  secondary  reparable  i 

li  mission  duration  during  which  secondary  reparable  i 
must  operate. 


If  perfect  knowledge  about  the  failure  rates  is  not  assumed 
and  the  Bayesian  approach  is  adopted  the  coefficients  are: 
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where, 

a,-  number  of  pacing  items  containing  secondary  reparable 
i 

Av  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (k^£a^) 

ni  number  of  observations  (data  points)  of  times  between 
failure  for  secondary  reparable  i 

’  total  observed  exposure  time  for  secondary  reparable  i 

[,  mission  duration  during  which  secondary  reparable  i 
must  operate 

empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(  /  ;  v.  [>  i  )  prior  distribution 

A 

Pi  empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(  {i  j_ )  prior  distribution. 

It  is  of  interest  to  note  that  for  situations  in  which 
the  airlift  model  is  applicable,  problem  P2  is  just  a 
special  case  of  the  generalized  assignment  problem.  In  the 
formulation  of  the  generalized  assignment  problem 
[Ref.  24:pp. 345-346] ,  each  agent  (or  man)  can  perform  or  be 
assigned  more  than  one  task  (or  job)  provided  that  the 
resource  available  to  the  agent  is  not  exceeded;  whereas,  in 
the  classical  assignment  problem  [Ref.  23:pp.  61-62],  each 

agent  can  be  assigned  to  only  one  task.  In  both 


formulations,  each  task  must  be  assigned  to  exactly  one  of 


the  agents.  The  airlift  model  can  be  interpreted  as  a 
special  case  of  this  generalized  assignment  problem  in  which 
one  agent  must  perform  all  N  tasks,  however,  he  has  ( +  1) 
alternative  ways  to  accomplish  each  task. 

Problem  P2  is  equivalent  to  PI  and  even  though  the 
number  of  variables  and  constraints  have  increased,  it  is 
the  preferred  formulation  because  it  eliminates  the  complex 
nonlinear  objective  function  of  PI.  The  coefficients  f j_j 
can  be  directly  computed  using  recursive  formulas  for  the 
Poisson  or  negative  binomial  cumulative  distribution 
functions;  whereas,  p^(x^)  could  not  since  it  was  itself  a 
function  of  the  decision  variable.  Additionally,  because  of 
the  assumption  of  the  independence  of  failures,  other 
failure  distributions  can  be  readily  introduced  into  this 
linear  formulation.  Furthermore,  the  formulation  is  general 
enough  so  that  additional  constraints  can  be  incorporated  as 
required  and  as  will  be  seen  in  the  following  chapter,  it 
lends  itself  easily  to  programming  in  an  algebraic  modeling 
language. 
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SOLUTION  METHODOLOGY  AND  IMPLEMENTATION 


The  purpose  of  this  study  is  to  determine  an  optimal 
stockage  level  of  spare  secondary  reparables  for  each  of  the 
three  different  scenarios  described  in  tho  introduction. 
Of  course  "optimality"  is  with  reference  to  the  failure 
model  adopted  and  the  relevance  of  the  data  used  to  estimate 
its  parameters.  The  determination  of  the  stockage  levels 
is  the  subject  of  this  chapter.  The  binary  linear  nature  of 
the  formulation  presented  in  the  preceding  chapter  suggests 
that  several  different  methods  of  solution  are  available  to 
determine  the  optimal  stockage  policy  for  a  given  scenario. 
Each  alternative  method  will  be  discussed  and  justification 
for  the  method  chosen  is  provided  in  addition  to  its 
implementation. 

A.  SOLUTION  APPROACHES  AND  METHODOLOGY 

The  mode  of  deployment  in  a  particular  scenario 
determines  the  number  of  constraints.  The  initial  solution 
approach  to  the  scenario  in  which  the  airlift  model  was 
applicable  utilized  the  technique  of  Lagrangian  relaxation. 
Although  it  is  standard  procedure  to  establish  bounds  for 
integer  programs  (such  as  problem  P2 )  using  a  linear 
programming  (LP)  relaxation,  it  was  initially  desired  to 
avoid  dependency  on  commercial  LP  solvers  not  locally 


,v:i 


available  to  the  combat  service  support  planners  of  the 
Second  Force  Service  Support  Group  at  Camp  Lejeune,  North 
Carolina.  Therefore,  Lagrangian  relaxation  was  adopted 
because  it  lends  itself  nicely  to  programming  in  readily 
available  FORTRAN  77.  It  must  be  emphasized  that  the 
structure  of  problem  P2  did  not  demand  that  this  specific 
solution  technique  be  applied.  In  fact,  the  LP  relaxation 
of  problem  P2  could  have  been  solved  using  an  LP  solver. 
Thus,  Lagrangian  relaxation  was  utilized  not  out  of 
necessity,  but  primarily  because  it  offered  a  realistic 
alternative  to  the  use  of  a  commercial  LP  optimizer. 

Briefly  stated,  the  technique  of  Lagrangian  relaxation, 
as  outlined  by  Fisher  [Ref.  25:pp.  10-21],  involves  moving 
complicating  constraints  into  the  objective  function  using 
the  product  of  the  Lagrangian  multiplier  and  the  constraint 
violation  as  a  penalty  term.  In  particular,  examination  of 
problem  P2  in  the  context  of  the  airlift  scenario  results  in 
the  following  binary  linear  formulation  (actually  a 
generalized  assignment  formulation): 

INDICES: 

/=1,...,.V  secondary  reparables  or  components 
j  —  0, ... ,  number  of  secondary  reparables  i 


DATA: 
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4  ln( p± ( j ) ) 

"V  total  shipping  weight  for  j  spares  of  type  i 

[  lbs} 

ACL  allowable  cargo  load  of  airplane  [lbs] 

ut  maximum  allowable  number  of  spares  for 

component  i 


DECISION  VARIABLE: 


-  f1 
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if  j  spares  of  type  i  are  selected 
otherwise 


FORMULATION: 
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The  only  complicating  constraint  is  the  allowable  cargo  load 
(ACL)  constraint.  The  relaxed  formulation  with  scalar 
Lagrangian  multiplier,  A  ,  consists  of  N  separable,  multiple 
choice  problems  which  are  easy  to  solve.  On  the  other  hand, 
for  those  situations  in  which  both  the  weight  and  volume 
constraints  of  problem  P2  must  be  considered  (i.e., 
Mediterranean  MAU  or  MPF  MAB  scenarios ) ,  the  Lagrangian 
multiplier,  X  ,  is  now  vector-valued  (i.e.,  A  =(  A]_,  A  2)  since 
two  constraints  must  be  relaxed,  and  this  is  no  longer  a 


trivial  problem  to  solve-  The  primary  reason  for  this 
complication  is  that  the  solution  approach  now  involves  a 
two-dimensional  search  which  is  more  difficult  to  program 
than  the  simple  bisection  search  associated  with  the  single, 
scalar  Lagrangian  multiplier.  The  problem  is  not  impossible 
to  solve  and  methods  exist,  such  as  a  coordinate  search 
method  suggested  by  DeWolfe,  Stevens,  and  Wood  [Ref. 26: 
pp.8-10],  which  appears  to  offer  good  results  for  this  type 
of  problem. 

As  additional  constraints  (such  as  budget,  holding 
(storage)  cost,  or  shortage  (stockout)  cost  constraints)  are 
added,  the  programming  difficulty  associated  with  the  search 
procedure  of  the  Lagrangian  function  becomes  prohibitive. 
Therefore,  to  accommodate  the  possible  inclusion  of 
additional  constraints  in  the  formulation  of  problem  P2,  it 
was  decided  to  utilize  the  Generalized  Algebraic  Modeling 
System  (GAMS)  in  conjunction  with  the  Zero/One  Optimization 
Methods  (ZOOM)  mixed  integer  program  solver.  The  use  of 
this  optimization  package  not  only  avoids  the  restrictive 
programming  considerations  discussed  above,  but  it  can  also 
incorporate  the  PEB  as  an  initial  step,  prior  to  solution  of 
P2  or  P3 . 

The  approach  to  the  PEB  is  to  maximize  the  marginal 
likelihood  function  associated  with  each  secondary  reparable 
i  with  respect  to  q  and  fi  i.  This  marginal  likelihood 


function  is  a  nonlinear  function  of  the  two  decision 


variables.  The  formulation  for  this  maximization  problem 
is : 


INDICES: 

,.V  secondary  reparable  or  component 


m  —  1 . ... ,  \ft  source  of  data  (e.g.,  tank  serial  number  508282) 


DATA: 


number  of  observations  of  component  i  failures 
for  data  source  m 


total  observed  exposure  time  for  component  i  of 
data  source  m. 


upper  bound  on  the  shape  parameter  of  the  Gamma 
prior  distribution  for  component  i. 


upper  bound  on  the  scale  parameter  of  the  Gamma 
prior  distribution  for  component  i. 


DECISION  VARIABLES: 


xi  shape  parameter  for  Gamma  prior  distribution 

associated  with  compenent  i. 


Pi  scale  parameter  for  Gamma  prior  distribution 

associated  with  component  i. 


FORMULATION: 
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Standard  maximization  by  taking  partial  derivatives  with 
respect  to  a ^  and  does  not  yield  a  satisfactory  closed 
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form  solution  which  can  be  practically  implemented.  Thus, 
this  nonlinear  maximization  problem  is  also  solved  using 
GAMS  in  conjunction  with  the  nonlinear  program  solver. 


MINOS  [Ref.  27].  These  values  are  then  incorporated  into 
the  optimization  of  problem  P2.  In  summary,  under  the 
current  solution  approach,  the  entire  model  can  be  coded  for 
computer  processing  on  either  a  mainframe  system  or  personal 
computer  and  solved  using  a  single  modeling  system  (i.e., 
GAMS)  and  its  compatible  solvers.  This  approach  is 
sufficiently  flexible  to  incorporate  additional  constraints 
or  decision  variables. 

B.  IMPLEMENTATION 

Since  the  decision  was  made  to  determine  the  optimal 
stockage  policy  for  all  scenarios  using  GAMS,  no  specifics 
of  the  original  solution  approach  involving  Lagrangian 
relaxation  will  be  mentioned  other  than  the  fact  that  it  is 
implemented  by  means  of  a  FORTRAN  77  computer  program.  The 
FORTRAN  program  is  based  upon  a  similar  one  developed  by 
DeWolfe  [Ref.  17:pp.  23-24].  A  copy  of  the  source  code 
developed  for  this  study  is  provided  in  Appendix  A  and 
examples  of  its  two  required  input  files  and  of  the  single 
output  file  are  given  in  Appendices  B  and  C,  respectively. 

The  solution  approach  adopted  centers  around  the  coding 
of  problems  P2  and  P4  for  computer  processing  using  the  GAMS 
language.  The  reader  is  referred  to  the  excellent  tutorial 
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by  Rosenthal  [Ref.  28]  for  the  specifics  of  this  versatile 
mathematical  modeling  language.  Upon  coding,  the  problems 
are  solved  using  the  GAMS  compatible  solvers  of  MINOS  and 
ZOOM  for  nonlinear  programs  and  mixed  integer  programs 
respectively.  Specifically,  the  source  code  provided  in 
Appendix  D  incorporates  the  successive  solve  feature  of  GAMS 
to  first  solve  the  nonlinear  maximization  problem 
represented  by  problem  P4.  The  results  of  this  maximization 

A  * 

provide  the  estimates  of  the  shape  (  )  and  scale  (  p ) 
parameters  of  the  prior  distributions  which  are  required 
before  the  second  binary  linear  maximization  ( P2 )  can  be 
performed.  The  final  result  is  the  optimal  stockage  level 
of  spare  secondary  reparables  for  the  given  scenario.  Even 
though  the  source  code  in  Appendix  D  is  specific  to  the 
airlift  model  ( P3 ) ,  GAMS  makes  it  very  easy  to  incorporate 
additional  constraints.  For  example,  if  a  volume  constraint 
is  also  applicable,  a  single  GAMS  PARAMETER  statement  to 
accomodate  the  volume  data  and  EQUATION  statement  to  reflect 
the  volume  constraint  are  the  only  additions  necessary. 

The  source  code  provided  is  compatible  with  either 
mainframe  or  personal  computer  usage.  When  run  on  the  IBM 
3033AP  mainframe  computer,  GAMS  version  2.05  was  utilized  in 


conjunction  with  MINOS  Version  5.1  and  ZOOM  version  2.1  and 
on  the  Zenith  Z-248  personal  computer  (equipped  with  80287-8 
math  coprocessor  and  640  Kbyte  memory  expansion  board),  PC 


GAMS  Version  2.04  was  utilized  in  conjunction  with  MINOS 
Version  5.0  and  ZOOM  Version  2.1. 


V.  RESULTS 


This  chapter  summarizes  results  obtained  using  actual 
operational  maintenance  data  to  determine  an  optimal  level 
of  spare  secondary  reparables  for  a  Marine  Amphibious  Unit 
(MAU)  in  a  realistic  exercise  scenario.  Verification  of  the 
model  is  accomplished  using  simulated  failure  data. 
Estimation  of  the  parameters  of  the  prior  distribution  by 
an  analysis  of  the  data  is  discussed  and  demonstrated;  this 
realizes  computationally  the  PEB  procedure.  Finally, 
computational  experience  for  both  the  mainframe  and  Zenith- 
248  personal  computer  is  briefly  discussed. 

A.  MODEL  TESTING 

The  data  used  in  the  development  and  testing  of  all 
models  was  provided  by  the  Operations  Section  of  the 
Supported  Activities  Supply  System  (SASSY)  Management  Unit 
in  coordination  with  the  Maintenance  Information  Systems 
Management  Office  of  the  Second  Force  Service  Support  Group 
at  Camp  LeJeune,  North  Carolina.  This  data  was  extracted 
from  the  ERO  History  File  of  the  Marine  Corps  Integrated 
Maintenance  Management  System  and  analyzed  according  to  the 
procedure  presented  in  Chapter  II.  This  analysis  provided 
the  basic  parameters  for  equation  (3.5)  (or  (3.4))  from 


which  the  objective  function  coefficients  were  generated 
using  GAMS . 

The  airlift  model  ( P3 )  was  used  to  determine  the  optimal 
stockage  policy  for  twenty-five  different  secondary 
reparables  associated  with  the  pacing  items  of  a  Marine 
Amphibious  Unit  participating  in  a  realistic  short-term 
exercise.  This  problem  generates  190  binary  variables  and 
26  constraints.  The  recommended  stockage  level  for  the 
secondary  reparables  associated  with  the  three  pacing  items 
(i.e.,  M60A1  Tank,  LTP7A1  Amtrack,  and  M101A1  105mm 
Howitzer)  in  this  exercise  is  presented  in  Table  2.  Table  2 
displays  the  levels  computed  using  both  the  GAMS/  MINOS/ZOOM 
package  and  the  Lagrangian  relaxation  technique. 

For  this  exercise  scenario  99.89%  of  the  allowable  cargo 
load  was  utilized  and  the  objective  function  value  was 
within  0.07%  of  the  upper  bound  established  during  the 
branch  and  bound  procedure  used  by  the  ZOOM  solver.  The 
probability  of  mission  completion  given  the  stockage  level 
was  0.4149  and  this  highlights  one  of  the  advantages 
associated  with  using  this  model.  This  low  probability 
implies  that  the  allowable  cargo  load  must  be  increased.  The 
only  way  to  accomplish  this  is  to  increase  the  number  of 
cargo  planes  assigned  to  this  exercise.  This  provides 
Marine  Corps  planners  with  a  tangible  basis  on  which  to 


TABLE  2 

MAU  AIRLIFT  SCENARIO  SPARE  STOCKAGE  LEVELS 
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justify  their  requests  for  additional  airlift  assets  from 
the  Military  Airlift  Command  of  the  U.S.  Air  Force.  Table  3 
summarizes  the  results  for  both  solution  approaches. 


TABLE  3 


SUMMARY  OF  RESULTS  FOR  MAU  AIRLIFT  SCENARIO 


Objective  Probability 
Solution  Function  of  Mission 
Techniaue  Value  ComDletion 


Deviation 

Percentage  from 

of  ACL  Used  Ontimalit 


GAMS /MINOS/ 

ZOOM  -0.8797 


0.4149 


99.89% 


0.07% 


Lagrangian 
Relaxation  -0.8820 
( w/Heuristic ) 


0.4139 


100% 


0.47% 


B.  MODEL  VERIFICATION 

Verification  of  the  modeling  approach  was  conducted  by 
using  the  model  to  verify  the  fact  that  the  negative 
binomial  predictive  function  reduces  to  the  Poisson-Gamma 
function  as  the  number  of  failures  observed  (and  accordingly 
the  total  exposure  time)  approaches  infinity.  Thus,  the 
Poisson  function  provides  the  limiting  case  or  best  possible 
stockage  level  attainable  assuming  that  the  failure  rates 
for  each  secondary  reparable  are  precisely  known. 

To  accomplish  this  verification,  the  MAU  airlift  model 
was  run  using  the  Poisson  objective  function  (3.4)  with 
fixed  failure  rates.  The  resulting  stockage  levels 
represent  the  best  stockage  policy  that  can  be  achieved. 
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Then,  using  those  failure  rates,  twenty-five  different  sets 
of  simulated  failure  data  were  generated  for  each  of  10,  25, 
50,  100,  200  and  500  failures.  These  simulated  failure 
times  were  randomly  generated  using  the  probability  integral 
transform  to  obtain  exponential  failure  times  from  uniform 
(0,1)  deviates  provided  by  the  LLRANDOMII  Random  Number 
Generation  Package  [Ref.  29] .  The  Bayesian  analysis  of  this 
data  was  performed  resulting  in  the  negative  binomial 
predictive  distribution  and  the  objective  function,  equation 
(3.5).  A  noninformative  Gamma  prior  distribution  with  ~x  =1 
and  p  =1  was  employed  in  this  process.  Stockage  levels  for 
this  method  are  displayed  in  Table  4.  Note  that  about  100 
failure  observations  are  needed  before  the  stockage  levels 
approach  those  computed  using  the  actual  failure  rates 
(i.e.,  the  optimal  stockage  policy  given  perfect 
information).  Figure  1  is  a  plot  of  the  interquartile  range 
( IQR )  of  the  probabiltiy  of  mission  completion  associated 
with  each  failure  sample  size.  This  graph  demonstrates  that 
as  the  number  of  observed  failures  increases,  the 
variability  (represented  by  the  IQR)  in  the  probability  of 
mission  completion  decreases. 

In  addition  to  the  above  verification,  the  fact  that  the 
model  and  solution  approach  responded  predictably  to  a  wide 
range  of  different  failure  observations  and  exposure  times 
indicates  that  the  model  and  approach  tend  to  be  robust. 
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TABLE  4 


MAU  AIRLIFT  SCENARIO  SPARE  STOCKAGE  LEVELS 
SPARE  STOCKAGE  LEVEL 

NEGATIVE  BINOMIAL 

NSN _ 10  25  50  100  200  500  POISSON 

2815001245387  3332  22 

2520002241867  4443  33 

2815002395819  2222  22 

2520001549632  2223  33 

2530000886657  1112  22 

2520001184942  1112  21 

2920010135802  1111  11 

1010010703803  0011  11 

1025019388105  4444  44 

1010010720397  1111  11 

1025010708993  1001  11 

2815004303480  5444  44 

2520003973384  5555  55 

2520001443385  3444  44 

5805014593410  5555  55 

2520008949532  2222  22 

2910011714636  1222  22 

2530011509757  2222  22 

5865012207848  3223  33 

4320012035660  3222  22 

2920002317276  1111  11 

2530000886650  1111  12 

1005011855059  4434  44 

1005011086434  5555  55 

1005011457709  5555  55 
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BAYESIAN  AND  EMPIRICAL  BAYES  ANALYSIS 


Bayesian  analysis  and  particularly  PEB  suggest 


themselves  when  the  failure  rates  of  different  units  (e.g.. 


tank  engines  on  different  tanks  or  amtrack  transmissions  on 


different  amtracks )  are  sufficiently  alike,  but  individual 


experience  is  small,  so  that  "borrowing  strength"  by  a 


suitable  pooling  process  can  improve  individual  estimates. 


In  this  study,  the  failure  rate  associated  with  each 


secondary  reparable  (  i )  was  assumed  to  follow  a  Gamma 


prior  distribution.  This  prior  distribution  represents  the 


probability  that  the  unknown  value  of  the  failure  rate  lies 


in  various  regions  of  the  parameter  space  [Ref.  20:p.  260], 


In  other  words,  incorporating  the  prior  distribution  enables 


a  tighter  bound  on  the  unknown  value  of  the  failure  rate  to 


be  obtained  than  would  be  possible  by  using  the  data 


pertaining  to  that  rate  alone.  The  primary  disadvantage  of 


the  ordinary  or  classical  Bayesian  approach  is  that  the 


prior  distribution's  parameters  (i.e.,  and  /J^)  are 


assumed  to  be  known.  This  deficiency  is  overcome  in  this 


study  by  employing  a  PEB  approach,  in  which  the  entire  data 


set  (  (nim .  )  associated  with  each  secondary 


reparable  i  is  used  to  estimate  the  parameters  of  the  prior 


distribution.  These  estimates  are  then  used  in  the  standard 


Bayesian  analysis  to  compute  the  posterior  and  predictive 


distributions.  The  results  are  unavoidably  influenced  by 
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the  particular  choice  of  prior  used  (here  a  Gamma). 
Sentitivity  tests  have  not  been  conducted,  but  indication  is 
that  the  procedure  is  rather  insensitive. 

A  specific  indication  of  this  insensitivity  is  seen  by 
comparing  the  airlift  model  stockage  levels  obtained  when 
minimal  prior  information  is  included  to  the  levels  obtained 
when  the  empirical  Bayes  estimates  of  the  prior  distribution 
are  included.  The  results  of  this  comparison  are  summarized 
in  Table  5. 

The  levels  in  the  column  labeled  "NONINFORMATIVE  PRIOR" 
were  obtained  by  running  the  airlift  model  with  the  randomly 
generated  failure  times  (described  in  the  previous  section) 
associated  with  the  number  of  failures  sample  size  of  ten. 
A  "diffuse"  or  noninformat ive  prior  distribution 
(i.e.,  ai=  =  was  used.  A  diffuse  prior  distribution  is 
"informationless  in  a  relative  sense  and  means  only  that  it 
is  diffuse  relative  to  the  sample  information" . [Ref .  30:p. 
198]  The  same  airlift  model  was  then  run  using  empirical 
Bayes  estimates  for  the  prior  parameters  and  these  levels 
appear  in  the  column  labeled  "INFORMATIVE  PRIOR".  These 
levels  are  only  slightly  closer  to  the  Poisson  levels  than 
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TABLE  5 


EMPIRICAL  BAYES  ESTIMATION 

STOCKAGE  LEVELS 
NEGATIVE  BINOMIAL 


NSN 


2815001245387 

2520002241867 

2815002395819 

2520001549632 

2530000886657 

2520001184942 

2920010135802 

1010010703803 

1025019388105 

1010010720397 

1025010708993 

2815004303480 

2520003973384 

2520001443385 

5805014593410 

2520008949532 

2910011714636 

2530011509757 

5865012207848 

4320012035660 

2920002317276 

2530000886650 

1005011855059 

1005011086434 

1005011457709 


NONINFORMATIVE 

PRIOR 


3 

4 
2 
2 
1 
1 
1 
0 

4 
1 
1 

5 
5 
3 
5 
2 
1 
2 
3 

3 
1 
1 

4 

5 
5 


INFORMATIVE 

PRIOR 


2 

4 

2 

2 

1 

1 

1 

0 

4 
1 
1 

5 
5 

4 

5 
2 
1 
2 
3 

3 
1 
1 

4 

5 
5 
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are  those  computed  when  minimal  prior  information  is 
considered. 

Even  though  these  stockage  levels  are  only  slightly 
closer  to  the  Poisson  levels  for  this  specific  instance,  by 
simply  using  the  available  data  for  10  failures  to  compute 
empirical  Bayes  estimates  of  the  prior  distribution 
parameters,  spare  stockage  levels  comparable  to  those 
obtained  in  the  previous  section  using  failure  sample  sizes 
of  between  50  to  100  can  be  achieved.  This  is  particularly 
significant  since  we  have  no  control  over  the  number  of 
failures  observed  for  each  secondary  reparable  or  their 
exposure  times  and  would  not  usually  observe  this  number  of 
failures  (i.e.,  50-100)  during  the  reporting  period  covered 
by  the  ERO  History  File.  Also,  since  operational 
reliability  data  may  not  be  available  for  the  exact  location 
where  the  MAGTF  is  deploying,  pooling  the  data  from  sources 
which  have  operated  in  similar  environments  has  the 
advantage  of  providing  more  realistic  failures  rates  than  if 
this  prior  information  were  not  considered. 

D.  COMPUTATIONAL  EXPERIENCE 

So  as  to  keep  the  implementation  and  solution  of  the 
models  developed  in  this  study  as  general  and  flexible  as 
possible  existing  software  packages  (i.e.,  STATGRAPHICS , 
GAMS,  etc.)  were  utilized.  The  goal  was  to  obtain  solutions 
in  a  relatively  user-friendly  personal  computer  environment 
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requiring  a  minimum  of  operator  interaction.  Therefore, 
computational  efficiency  in  terms  of  computer  run  time  (or 
algorithmic  complexity)  was  not  of  primary  concern  and  is 
included  here  solely  for  completeness.  No  attempt  was  made 
to  employ  special  algorithms  that  exploit  problem  structure. 
This  is  not  to  downplay  the  importance  of  this  factor,  but 
rather  to  emphasize  that  during  this  study  the  primary 
concern  was  to  obtain  solutions  to  real-world  problems  as 
easily  and  efficiently  as  possible  with  a  minimum  of 
programming  effort. 

Running  time  for  the  airlift  model  on  the  IBM  3033AP 
using  the  GAMS/MINOS/ZOOM  package  was  3.9  seconds  in 
contrast  to  the  176  seconds  obtained  when  run  on  the  Zenith 
Z-248  PC.  It  is  possible  on  the  IBM  3033AP  to  combine  the 
PEB  calculations  with  the  maximization  of  problem  P3  by 
using  the  successive  solve  feature  of  GAMS;  however,  this 
cannot  be  done  on  the  Zenith  Z-248  PC  since  internal 
parameters  in  the  PC  version  of  GAMS  are  exceeded. 
Therefore,  each  of  the  empirical  Bayes  maximizations  must  be 
performed  separately  when  using  this  PC.  The  resulting 
estimates  have  to  be  manually  input  into  the  appropriate 
GAMS  PARAMETER  statements  before  maximization  of  problem  P3 
can  be  accomplished. 

Even  though  the  Lagrangian  relaxation  procedure  with 
heuristic  was  not  adopted  as  the  solution  technique  for  the 


reasons  discussed  in  Chapter  IV,  it  provides  much  faster 
solution  times  to  the  airlift  model  ( P3 )  on  both  the  IBM 
3033AP  and  Zenith  Z-248  PC.  The  Lagrangian  procedure,  coded 
in  FORTRAN  77,  provides  a  running  time  of  0.32  seconds  on 
the  IBM  3033AP  and  about  30  seconds  on  the  Zenith  Z-248  PC 
where  it  was  coded  using  Microsoft  FORTRAN  (Version  3.2). 

Finally,  comparison  of  the  stockage  levels  obtained 
using  GAMS/MINOS/ZOOM  with  those  obtained  using  LP 
relaxation  for  each  of  the  150  simulated  failure  time  cases 
revealed  that  in  every  case  the  LP  relaxation  provided 
identical  stockage  levels  with  the  exception  that  it 
fractionated  for  exactly  one  of  the  twenty-five  secondary 
reparables.  So,  if  a  general  policy  was  adopted  to  round 
down  to  the  lower  of  the  fractionated  stockage  levels  ( so  as 
not  to  violate  the  ACL  contrast),  the  LP  relaxation  would 
provide  comparable  results  to  those  obtained  using  the  mixed 
integer  program  solver,  ZOOM.  This  is  an  important  result 
because  the  run  times  associated  with  the  LP  relaxation  were 
significantly  less  than  those  using  ZOOM.  In  particular, 
if  the  rounding  down  policy  is  applied  to  the  LP  relaxation 
of  the  airlift  model  in  this  study,  the  resulting  integer 
stockage  level  provides  an  objective  function  value  which  is 
within  1.87%  of  the  LP  upper  bound  or  1.73%  of  the  upper 
bound  established  by  the  branch  and  bound  procedure  utilized 
by  ZOOM.  This  is  not  quite  as  good  as  the  result  obtained 


directly  from  ZOOM  (within  0.07%  of  optimality),  however, 
for  those  situations  in  which  computer  run  time  may  be  a 
consideration  (e.g.,  as  problem  size  increases)  or  in  which 
very  tight  bounds  may  not  be  required,  a  very  good  solution 
can  be  obtained  from  the  much  simpler  and  less  costly  LP 
relaxation. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  thesis,  the  problem  of  determining  spare 
stockage  levels  for  deployable  MAGTFs  was 
solved  by 

(a)  stochastically  modeling  demand  utilizing  operational 
data  from  the  Marine  Corps  Integrated  Maintenance 
Management  System,  and  then 

(b)  using  mathematical  programming  to  optimize  a 
suitable,  operationally  relevant  objective  function. 

To  accomplish  this,  a  new  measure  of  effectiveness  was 

defined  which  provides  the  operational  commander  on  the 

ground  with  an  objective  indicator  of  the  overall  impact  of 

the  stockage  policy  on  the  probability  of  mission 

completion.  This  measure  of  effectiveness  was  evaluated  as 

the  end  result  of  a  probabilistic  modeling  step,  a 

statistical  data  analysis,  and  a  subsequent  optimization. 

The  generality  of  the  model  formulation  makes  it 
adaptable  to  any  of  the  three  scenarios  of  present  interest 
to  the  Marine  Corps.  The  use  of  a  PC-compatible  algebraic 
modeling  language  (GAMS)  makes  it  particularly  easy  to  add 
or  delete  constraints  as  required.  Specifically,  the  model 
and  solution  approach  were  applied  to  a  realistic  exercise 
scenario  and  an  essentially  optimal  stockage  policy  was 
determined  which  maximized  the  probability  of  mission 


completion  within  the  constraints  imposed  by  the  deployment 
mode . 

During  the  analysis  of  the  MIMMS  data  base  ( i . e . , ERO 
History  File)  numerous  problems  and  inconsistencies  were 
encountered.  As  a  result,  the  estimates  of  the  failure 
rates  for  some  of  the  secondary  reparables  are  higher  than 
would  be  the  case  if  all  the  potential  times  to  failure 
could  be  extracted  from  the  data  base.  To  overcome  this 
discrepancy,  it  is  highly  recommended  that  the  quarterly 
update  procedure  described  in  Chapter  II  be  conducted  at  the 
Force  Service  Support  Group  level.  If  this  cannot  be  done, 
then  at  the  very  least  the  meter  readings  associated  with 
all  pacing  items  should  be  recorded  and  input  into  MIMMS 
prior  to  each  quarterly  update.  This  would  insure  that 
beginning  and  ending  meter  readings  are  available  for  each 
serialized  pacing  item  during  the  reporting  period  covered 
by  the  ERO  History  File  and  this  would  prevent  the  loss  of 
valuable  exposure  time  data. 

The  limited  strategic  lift  assets  available  to  deploying 
Marine  Corps  units  demands  that  effective  utilization  be 
made  of  these  resources.  Since  only  a  small  percentage  of 
these  assets  can  be  allocated  for  spare  parts  inventories, 
emphasis  must  be  placed  on  stocking  those  parts  which  would 
have  the  greatest  impact  on  the  completion  of  the  assigned 
mission.  The  reliability-based  model  developed  in  this 
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LISTING  OF  SOURCE  CODE  FOR  LANGRANGIAN  RELAXATION 


PROGRAM  RELAX2 

CHARACTER  NSNXl 3 , N0MEN*7 

INTEGER  NCOMP, NEQUIP, NSPARE, MAXSPR, IDNUM, BSPARE, SPBND, NDENS, 
XNOPER , NXPOSE 

REAL  *8  P,W, MAXWT, DINF, OPTDEV , ZLOW, ZUP , PSUC, HI LAMD, WT , MTBF, 
XDURMIS, ZHLOW, XPOSE, ALPHA, BETA 

PARAMETER  C NCOMP =25 , NEQUIP=3 , MAXSPR=1 0 , MAXWT=44400 .  0D0, 

XDINF=1 . 0D20) 

DIMENSION  NSN (NCOMP), NOMEN ( NCOMP ), NSPARE( NCOMP ) ,IDNUM( NCOMP), 
xP ( NCOMP, 0= MAXSPR) , W( NCOMP , 0 : MAXSPR ) , BSPAREI NCOMP ) , SPBND( NCOMP), 
X, NDENS ( NCOMP ) , NOP ER( NCOMP ), XPOSEC NCOMP ), NXPOSEC NCOMP ) , 

*DURMIS< NCOMP), ALPHA (NCOMP), BETA (NCOMP) 

CALL  INIT2(DINF,MAXSPR, N EQUIP , MAXWT , P , W , NSN , NOMEN , HI L AMD , NSPARE, 
XSPBND,  NCOMP, IDNUM, NDENS, NOPER,XPOSE, NXPOSE, DURMIS, ALPHA, BETA) 

CALL  30UND(DINF,P,W, NCOMP, MAXSPR, SPBND, MAXWT, HI L AMD, NSPARE, ZLOW, 
*ZUP,WT, BSPARE) 

CALL  HERIST(DINF,P,W, NCOMP , MAXSPR , ZH LOW, WT, MAXWT , SPBND, DSPARE, 
XLARRAY) 

WRITE  (10,50) 


50  FORMAT  ( '1 '/'O' , • - 

*« - »  ) 

OPTDEV  =  -1.0D0  x  100. 0D0  *  (1.0D0  -  (ZHLOW/ZUP)) 


WRITE  (10,100)  ZHLOW 

100  FORMAT  ( IX, 'OBJECTIVE  FUNCTION  VALUE  IS  :',D17.10) 

WRITE  (10,150)  OPTDEV 

150  FORMAT  ('O', 'THIS  SOLUTION  IS  WITHIN  ',D17.10,'  X  OF  THE  OPTIMAL') 
PRINT  x, *ZUP=' ,ZUP, 'ZLOW= ' ,ZLOW, ' ZHLOW= ' , ZHLOW 
P3UC  =  EXP (ZHLOW) 

WRITE  (10,175)  PSUC 

175  FORMAT  ( 'O' , 'P(MISSION  SUCCESS)  =  ’,F7.5) 

WRITE  (10,180) 

130  FORMAT  (  '0'  ,  ' - ', 

x  ' - «  ) 

WRITE  (10,200) 

200  FORMAT  ( ' 0 ’ , 5X, ' NSN ' , 1  OX , ' END  ITEM' , 6X, *  SPARE  LEVEL') 

DO  250  1=1, NCOMP 

WRITE  (10,300)  NSN( I ) , NOMEN ( IDNUM( I ) ) , BSPARE( I ) 

250  CONTINUE 

300  FORMAT  ( • 0 ’ , A13, 7X, A7 , 7X, 14) 

PRINT  x, 'WEIGHT  USED= ', WT, ' MAXWT= ', MAXWT 

STOP 

END 


SUBROUTINE  INIT2C DINF, MAXSPR, NEQUIP , MAXWT , P , W, NSN , NOMEN , HILAMD, 
*NSPARE,SPBND,NCOMF, IDNUM, NDENS, NOPER, XPOSE, NXPGSE, DURMIS, ALPHA, 
XBETA) 

INTEGER  NCOMP 

REAL  *8  MAXWT 

INTEGER  I DNUMC  NCOMP ) , NDENSC NCOMP ) , NOPER ( NCOMP ) , LMT , MAXSPR, 
*NSPAREC NCOMP ) , SPEND (NCOMP ), NEQUIP, NXPOSEC NCOMP) 

REALMS  WC NCOMP, 0: MAXSPR), DURMIS ( NCOMP ) , DINF, PI , Q , BETAC NCOMP ) , 
*P TEMPI , PC NCOMP, 0 : MAXSPR) , HILAMD, HLTEMP, CUMSUM, XPQS EC NCOMP) , PP, 
XALPHACNCOMP) 

CHARACTER  NSNX13C NCOMP ), N0MENX7 ( NCOMP) 

READ  (1,10)  (NSN(I),IDNUM(I),XPOSE(I), NXPOSEC I),W(I,1),SPBND(I), 
XALPHAC I ) , BETAC I ) ,  1=1, NCOMP) 

10  FORMAT  CA13,2X,I2,2X,F9.1,2X,I4,2X,F6.1,2X,I3,2X,F7.4,2X,F8.2) 

READ  C  2 , 20 )  ( NOMEN ( I ) , NDENS C I ) , NOPERC I ) , DURMIS C I ) ,  1  =  1, NEQUIP) 

20  FORMAT  ( A7 , 4X, 13, 4X, 13, 4X, F6  . 1 ) 

DO  21  1=1, NCOMP 

PRINT  X. NSN Cl), I DNUMC I), XPOSECI) , NXPOSEC I),W(I,  1),SPBNDCI), 
XALPHACI).BETA(I) 

21  CONTINUE 

DO  22  1=1, NEQUIP 

PRI N'T  x ,  NOMEN ( I ) ,  NDENSC  I ),  NOPERC  I ) ,  DURMISC I ) 

22  CONTINUE 


DO  25  1=1, NCOMP 
NSPAREC I )  =  0 
DO  5C  J=0, MAXSPR 
PCI,J)  =  DINF 

IF  CJ  .NE.  1)  WCI,J)  =  DINF 
5C  CONTINUE 

25  CONTINUE 

DO  75  1=1, NCOMP 

PI  =  CXPOSECD  +  BETACI) )/( NDENSC IDNUMCI )) * 

XDURMISCIDNUMCI) )  +  XPOSECI)  +  BETACI ) ) 

Q  =  1 . 0D0  -  PI 

PP  =  PI  xx  ( NXPOSEC I )  +  ALPHA(I)) 

LMT  =  NDENSC I DNUMC I ) )  -  NOPERC I DNUMC I ) ) 

DC  100  J  =  0 , 5PBNDC I ) 

IF  CJ  .NE.  1)  W(I,J)  =  J*W(I,1) 

CUMSUM  =  0.0D0 
DO  125  K=  0 , ( LMT  +  J) 

IF  CK  .EQ.  0)  THEN 
PTEMP1  =  1.0D0 
ELSE 

PTEMP1  =  PTEMP1X(CCCNXPOSE(I)+ALPHACI))+(K-1))X1.0DO)/K)XQ 
END  IF 

CUMSUM  =  CUMSUM  +  PTEMP1 
125  CONTINUE 

CUMSUM  =  CUMSUM  *  PP 
P  C I , J )  =  LOGC  CUMSUM) 

PRINT  X, *I=',I, ’J=',J, ,LN=',P(I,J), 'CUMSUM=' , CUMSUM 
100  CONTINUE 

75  CONTINUE 


CJ  (M 


HI L AMD  =  -DINF 
DO  200  1=1,  NCOMP 

DO  225  J=1,SPBND(I) 

HLTEMP  =  (  P(I,J)  -  PCI,0)  )  /  WCI,J) 

IF  (HLTEMP  .GT.  HILAMD)  HI L AMD  =  HLTEMP 
25  CONTINUE 
CO  CONTINUE 

PRINT  x, 1 HILAMBDA= 1 , HILAMD 

RETURN 

END 


SUBROUTINE  BOUND( DINF, P , W, NCOMP , MAXSPR , SPBND, MAXWT , HI L AMD, NSPARE 
XZLOW, ZUP , WT , BSP ARE ) 


INTEGER  NCOMP, BSPAREC NCOMP), MAXSPR, SPBNDC NCOMP ),NSPARE( NCOMP) 


REAL *8  BSTL AM, ZLOW, WT, P( NCOMP, 0 : MAXSPR ) , WC NCOMP, 0 : MAXSPR) , DINF 

*MAXWT, HILAMD, EPS, LFTEND, RTEND, ZLBEST , LAMBDA, 


EPS  =  HILAMD/1 . 0D04 
LFTEND  =  O.ODO 
RTEND  =  1 . 1D0XHILAMD 
ZLBEST  =  -DINF 

50  LAMBDA  =  (LFTEND  +  RTEND)/2.0D0 

CALL  MAXFCK DINF, P,W, NCOMP, MAXSPR, SPBND, MAXWT, LAMBDA, NSPARE, ZLOW 
*ZUP,WT) 

IF  (WT  .LE.  MAXWT)  THEN 
RTEND  =  LAMBDA 
IF  (ZLOW  .GE.  ZLBEST)  THEN 
ZLBEST  =  ZLOW 
BSTL AM  =  LAM3DA 
END  IF 
ELSE 

LFTEND  =  LAMBDA 
END  IF 

IF  (  (RTEND  -  LFTEND)  . GT .  EPS  )  GOTO  50 

CALL  MAXFEA( DINF, P , W, NCOMP, MAXSPR, SPBND, BSTLAM,BSPARE,WT, ZLOW) 
RETURN 


SUBROUTINE  MAXFCT ( DINF, P , W, NCOMP , MAXSPR, SPBND, MAXWT , LAMBDA, NSPARE 
*ZLOW,ZUP,WT) 

INTEGER  NCOMP, NSPAREC NCOMP), MAXSPR, SPBND( NCOMP), INDEX 


REAL  *8  LAMBDA, ZLOW.WT, PC NCOMP, 0 : MAXSPR ) ,WC NCOMP, 0 : MAXSPR ) , DINF, 
XCTOTAL , CMPMAX, ZVALJ, OBJECT, OBJVAL,ZUP, MAXWT 


WT  =  O.ODO 
CTOTAL  =  O.ODO 
ZLOW  =  O.ODO 


,%  - 


50 


25 


DO  25  1=1, NCOMP 
CMPMAX  =  -DINF 
INDEX  =  0 

DO  50  J=0, SPBKDC I ) 

ZVALJ  =  P(I,J)  -  LAMBDAXWCI, J) 
OBJ VAL  =  P(I,J) 

IF  (ZVALJ  .GE.  CMPMAX)  THEN 
NSPARE(I)  =  J 
CMPMAX  =  ZVALJ 
03JFCT  =  OBJ VAL 
INDEX  =  J 
END  IF 
CONTINUE 

CTOTAL  =  CTOTAL  +  CMPMAX 
ZLOW  =  ZLOW  +  OBJFCT 
WT  =  WT  +  W( I, INDEX) 

CONTINUE 


ZUP  =  CTOTAL  +  LAMBDA*MAXWT 

RETURN 

END 

SUBROUTINE  MAXFEA ( DINF, P, W, NCOMP, MAXSPR, SPBND, BSTL AM, BSPARE, WT , 
XZLOW) 

INTEGER  NCOMP, BSPAREC NCOMP ), MAXSPR, SPBND( NCOMP) , INDEX 

REAL  *8  BSTL AM, ZLOW.WT, PC NCOMP, 0 : MAXSPR ), W( NCOMP , 0 : MAXSPR) , DINF, 
KCMPMAX, ZVALJ, OBJFCT, OBJ VAL 


WT  =  O.ODO 
ZLOW  =  O.ODO 

DO  25  1  =  1,  NCOMP 
CMPMAX  =  -DINF 
INDEX  =  0 


tv 


50 

25 


DO  50  J  =  0 , SPBNDC I ) 

ZVALJ  =  P(I,J)  -  BSTLAM*W(I, J) 
OBJ VAL  =  P C I ,  J  ) 

IF  (ZVALJ  .GE.  CMPMAX)  THEN 
BSPAREC I )  =  J 
CMPMAX  =  ZVALJ 
OBJFCT  =  OBJ VAL 
INDEX  =  J 
END  IF 
CONTINUE 

ZLOW  =  ZLOW  +  OBJFCT 
WT  =  WT  +  W( I, INDEX) 

CONTINUE 


RETURN 

END 


SUBROUTINE  HERISTC DINF, P , W , NCOMP , MAXSPR, ZHLOW, WT , MAXWT , SPBND, 
XBSPARE, L ARRAY ) 

INTEGER  NCOMP , BS PARE ( NCOMP ), MAXSPR, SPBNDC  NCOMP) , INDEX 

REAL *8  DINF, P( NCOMP, 0 : MAXSPR ) ,W( NCOMP, 0 : MAXSPR ) ,WT1,WT, MAXWT, 
XXSSWT , CHKVAL , L ARRAY (NCOMP ) , NUMER, DENOM, L BEST  ,  ZHLOW 

CHKVAL  =  DINF/1.1D0 
HT1  =  WT 

XSSWT  =  MAXWT  -  WT1 


(NUMER 


NUMER/ DENOM 


0.0D0)  LARRAYC I )  =  -L ARRAY! I ) 


DO  50  1=1, NCOMP 

IF  ( ( BSPAREC I )  .GE.  SPBNDCI))  .OR.  C  PC  I , CBS 
XCHKVAL))  THEN 

L ARRAY ( I )  =  -DINF 
GOTO  50 
END  IF 

NUMER  =  P(I,(BSPARE(I)+1)>  -  PC  I , BSPARECI ) ) 
DENOM  =  W( I , (BSPARECI )+l) )  -  W( I , BSPAREC I ) ) 
IF  ((DENOM  .GT.  XSSWT)  .OR.  (NUMER  .LE.  1.0 
L ARRAY ( I )  =  -DINF 
GOTO  50 
ELSE 

LARRAYC I)  =  NUMER/DENOM 
END  IF 

IF  ( L ARRAY ( I )  , LT .  O.ODO)  LARRAYC I )  =  -LARR 
CONTINUE 
LBEST  =  -DINF 
DO  75  1=1, NCOMP 

IF  ( L ARRAY ( I )  . GT .  LBEST)  THEN 
LBEST  =  LARRAYC I ) 

INDEX  =  I 
END  IF 
CONTINUE 

IF  (LBEST  .LT.  (-CHKVAL))  GOTO  100 

WT1  =  WT1  -  W( INDEX, BSPAREC INDEX) )  +  W(INDEX,( 

BSPAREC INDEX)  =  BSPAREC INDEX)  +  1 

IF  ( WT1  .LE.  MAXWT)  GOTO  25 

ZHLOW  =  O.ODO 

WT  =  O.ODO 

DO  125  1=1, NCOMP 

ZHLOW  =  ZHLOW  +  P ( I , BSPAREC I ) ) 

WT  =  WT  +  WC I, BSPAREC I)  ) 

CONTINUE 

RETURN 

END 


(PCI,  CBSPARECD  +  l)) 


1 . QD-1 0 ) ) 


;i)  .GT.  LBEST)  THEN 
LARRAYC I ) 

I 


WC INDEX, (BSPAREC INDEX)+1) ) 


rv  rwnn  nji  n  w  •ui  njr  —  ■  n 


APPENDIX  B 
SAMPLE  INPUT  FILES 


PRINCIPAL  END  ITEM  FILE 


M60A1  5  4  60.0 
LVTP7A1  12  11  40.0 
M101A1  5  4  2000.0 


TABLE  B-l 

PRINCIPAL  END  ITEM  FILE  SPECIFICATIONS 


COLUMN 

DATA  DESCRIPTION 

DATA  TYPE 

1-7 

Nomenclature  of  end  item 

Character 

8-11 

Blank 

12-14 

Quantity  of  end  item 

Integer 

15-18 

Blank 

19-21 

Number  of  end  items 
required  operational 

Integer 

22-25 

Blank 

26-31 


Mission  duration 


Real*8 


SECONDARY  REPARABLE  COMPONENT  FILE 


2815001245387 

1 

2091.0 

7 

2520002241867 

1 

2091 . 0 

8 

2315002395819 

1 

2091.0 

3 

2520001549652 

1 

2091  .  0 

5 

2530000836657 

1 

2091.0 

1 

2520001184942 

1 

2091  ,  0 

1 

2920010135802 

1 

2091  .  0 

1 

10100] 0703303 

1 

2091 . 0 

1 

1025019388105 

1 

2091.0 

5 

1010010720397 

1 

2091.0 

3 

1025010708993 

1 

2091.0 

1 

2815004303480 

o 

3253 . 0 

8 

2520003973384 

2 

3258.0 

23 

2520001443385 

2 

3258 . 0 

14 

5305014593410 

2 

3258 . 0 

10 

2520008949532 

2 

3258 . 0 

3 

2910011714636 

2 

3258.0 

1 

2530011509757 

2 

3258 . 0 

4 

5865012207848 

2 

3253 . 0 

2 

4320012035660 

2 

3258.0 

3 

2920002317276 

2 

3258 . 0 

1 

2530000886650 

2 

3258 . 0 

1 

1005011355059 

3 

20996 . 0 

3 

1005011086434 

3 

20996 . 0 

11 

1005011457709 

3 

20996 . 0 

7 

2000 . 0 

5 

0.34000 

58.6903 

1000.0 

5 

1  .  00000 

1 .0000 

500 . 0 

5 

1  .  00000 

1 . 0000 

500 . 0 

5 

0.10110 

13.1850 

100.0 

10 

1  .  00000 

1 .0000 

250 . 0 

5 

1.00000 

1.0000 

500 . 0 

10 

1.00000 

1.0000 

1000 , 0 

10 

1  .  00000 

1  .0000 

100.0 

5 

1  .  00000 

1  .  0000 

2000 . 0 

5 

1.00000 

1.0000 

1  000.0 

5 

1  .  UOOOO 

1.0000 

500 . 0 

5 

0 . 44330 

142 .9404 

10QQ  .  0 

5 

0  .69700 

69 . 0598 

2000 , 0 

5 

1 . 00000 

1.0000 

250 . 0 

5 

0 . 20310 

27.7343 

500 . 0 

5 

1  .  00000 

1  .0000 

1QQ  .  0 

10 

1  .00000 

1 .0000 

1000  .  0 

10 

1 . 00000 

1  .0000 

100.0 

10 

1 . 00000 

1  .  0000 

500 . 0 

10 

1  .00000 

1  .0000 

500 . 0 

10 

1  .00000 

1.0000 

250.0 

5 

1 . 00000 

1.0000 

1000.0 

5 

1 . 00000 

1  .  0000 

250.0 

5 

1  .00000 

1  .0000 

500.0 

5 

1.00000 

1.0000 

TABLE  B-2 


SECONDARY  REPARABLE/COMPONENT  FILE  SPECIFICATIONS 

COLUMN  DATA  DESCRIPTION  DATA  TYPE 


1-13 

National  Stock  Number  (NSN) 

Character 

14-15 

Blank 

16-17 

Identification  Number 
(1  for  M60A1,  2  for  LVTP7A1 , 

3  for  M101A1 ) 

Integer 

18-19 

Blank 

20-28 

Total  exposure  time  for 
component 

Real*8 

29-30 

Blank 

31-34 

Number  of  observed  failures 
for  component 

Integer 

35-36 

Blank 

37-42 

Component  shipping  weight 

Real*8 

43-44 

Blank 

45-47 

Maximum  possible  stockage 
for  component 

Integer 

48-49 

Blank 

50-56 

Shape  parameter  of  prior 
distribution  for  component 
failure  rate 

Real*8 

57-58 

Blank 

59-66 

Scale  parameter  of  prior 
distribution  for  component 
failure  rate 

Real*8 

APPENDIX  C 


SAMPLE  OUTPUT  FILE 

OBJECTIVE  FUNCTION  VALUE  IS  : -0 . 8S20263158D+Q0 

THIS  SOLUTION  IS  WITHIN  0 . 476 1377 367 D+0 0  X  OF  THE  OPTIMAL 

PCMISSION  SUCCESS)  =  0.41394 


IT 

NSN 

END  ITEM 

SPARE  LEVEL 

2815001245387 

M60A1 

2 

u. 

2520002241867 

M6  0A1 

3 

—r 

2S15002395819 

M60A1 

2 

La 

2520001549632 

M60A1 

2 

z 

r* 

2530000386657 

M60A1 

3 

Ui 

> 

0 

2520001184942 

M60A1 

2 

0 

2920010135802 

M60A1 

1 

H 

< 

101C0107  03803 

M60A1 

1 

Q 

LU 

U 

1025019388105 

M60A1 

5 

3 

3 

1C1001Q720397 

M60A1 

1 

0 

ft 

a. 

1025010708993 

M60A1 

1 

lii 

q: 

2815004303480 

LVTP7A1 

4 

2520003973334 

LVTP7A1 

5 

2520001443385 

LVTP7A1 

4 

5805014593410 

LVTP7A1 

5 

2520008949532 

LVTP7A1 

2 

2910011714636 

LVTP7A1 

3 

2530011509757 

LVTP7A1 

2 

5365012207348 

LVTP7AI 

3 

4320012035660 

LVTP7A1 

2920002317276 

LVTP7A1 

1 

2530000886650 

LVTP7A1 

2 

1005011855059 

M101A1 

4 

1005011086434 

M101A1 

5 

1005011457709 

M101A1 

5 

78 


/ 


APPENDIX  D 


WM 

IvLsI 


S 


r 


«• 

$ 


9 


80 


GAMS  SOURCE  CODE 


end  item 
component 


K  spare  stockage  level  /S0*S1Q/ 


/E1875,E0346,E0640/ 
/ 1*25/ 


dummy  index 


/D0*D11/; 


SCALARS 

WTCAP  total  weight  capacity  (lbs)  /4440Q/ 
PSUC  probability  of  mission  completion; 


** 

PARAMETERS 

u. 

EQDENS(I) 

density 

/E1875 

5 

E0S46 

12 

u. 

E0640 

5/; 

<1 

2 

PARAMETERS 

2 

NOPER(I) 

number  er 

u J 
> 

/E1375 

4 

o 

E0846 

11 

0 

E0640 

4/; 

h- 

<( 

PARAMETERS 

DURMIS ( I ) 

mission 

0 

Ui 

u 

/E1875 

60 

E0846 

40 

a 

E0640 

2000/; 

a 

0 

PARAMETERS 

£t 

ft 

SHIPWT ( J ) 

shipping 

uj 

/I 

2000 

fc 

2 

1000 

3 

500 

4 

500 

5 

100 

6 

250 

7 

500 

8 

1000 

9 

100 

10 

2000 

11 

1000 

12 

500 

13 

1000 

14 

2000 

15 

250 

16 

500 

17 

100 

18 

1000 

19 

100 

20 

500 

21 

500 

22 

250 

23 

1000 

24 

250 

25 

500/; 

required  operational  during  mission 


duration  for  end  item  i  (hours  rounds  miles) 


t-tts 


5 


3 

Ui 

u 

G 

Q 

O 

E 

o. 

lit 

C 


PARAMETERS 

SPRMAX(J) 

/I 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 
17 
13 

19 

20 
21 
22 

23 

24 

25 

PARAMETERS 
NOBSC J) 

/I 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


max  number  of  spares  allowed  for  component  3 

5 

5 

5 

10 

5 

10 

10 

5 

5 

5 

5 

5 

5 

5 

5 

10 

10 

10 

10 

10 

5 

5 

5 

5/; 


number  of  data  points  observed  for  component  j 

8 

3 

5 

1 

1 

1 

1 

5 

3 

1 

8 

23 

14 

10 

3 
1 

4 
2 
3 
1 
1 
3 

11 

7/; 


lV 


kyfil 

1 


»5 


f 


& 

i 

;s 

$ 

*!£»! 


i 

Is 


PARAMETERS 
ALPHA ( J ) 
/I 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


PARAMETERS 
BETAC J) 
/I 
2 

3 

4 

5 

6 
7 


shape  paramtr  for  Gamma  prior  distribution 
0,34 
1 
1 

0  .1011 
1 
1 
1 
1 
1 
1 
1 

0.4433 

0.6970 

1 

0.2031 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1/; 


scale  paramtr  for  Gamma  prior  distribution 
58.6903 
1 
1 

13.1850 

1 

1 

1 


142.9404 
69 . 0598 
1 

27.7348 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1/; 


PARAMETERS 

OK(I.J)  equals  1  if  component  j  is  a  component  of  end  item  i 
/E1875 . ( 1*11 )  1 

E0846 .( 12X22)  1 

£0640.(23X25)  1/; 

PARAMETERS 

TXPOSE(I)  total  observed  exposure  for  end  item  i  (hrs  rds  mi) 
/E1875  2091 

E0S46  3258 

E0640  20996/; 

PARAMETER 

LMT(I,J)  max  number  of  failures  for  component  j; 

LMT ( I , J ) $ (  0K(I,J)  EQ  1  )  =  EQDEHS(I)  -  NOPER(I); 

SETS 

OKLVL 1 ( I ,  J ,  K )  feasible  spare  stockage  levels  for  component  j; 
0KLVL1(I,J,K)$((GRD(K)-1  LE  SPRMAX(J))  AND  (0K(I,J)  EQ  1))  =  YES 

x  The  following  set  is  necessary  to  obtain  the  cumulative  distributi 
x  function  of  the  negative  binomial  density. 


0KLVL2(I, J,K,l); 

CKLVL2(I, J,K,L)$(  (ORD(K)-l  LE  SPRMAX(J))  AND  (OK(I.J)EQ  1)  AND 
(ORD(L)-l  LE  LMT(I,J)+0RD(K)-1)  )  =  YES; 

DISPLAY  0KLVL2; 

OPTIONS  DECIMALS  =  8; 

PARAMETER 

P(I*J)  negative  binomial  paramtr 

Q(I,J)  one  minus  the  p  paramtr 

PP(I.J)  p  paramtr  raised  to  power  of  nobs  plus  alpha 

TWTCI.J.K)  total  snipping  wt  of  k  spares  of  component  j 

C(I,J,K)  cumulative  density  for  negative  binomial 

F(I,J,K)  natural  logarithm  of  cumulative  density 

PTMPC I, J , K, L )  temporary  paramtr  needed  to  sum  cdf  recursively; 

P( I »  J ) $( 0K( I >  J  )  EQ  1)  = 

(TXPOSE(I)  +  BETA(J))/(( EQDENSC I ) XDURMISC I ) )  +  TXPOSECI)  +  BETA(J) ) 
Q(I,J)$(0K(I,J)  EQ  1)  =  1  -  P( I » J ) ; 

PP(I, J)$(0K(I, J)  EQ  1)  =  P(I,J>  x*  (NOBS(J)  +  ALPHAC  J ) ) ; 

TWT(I, J,K)$0KLVL1(I, J,K)  =  SHIPWT(J)  X  (ORD(K)  -  1); 

PTMP(I, J,K,"DO")$OKLVL1(I, J,K)  =  1; 

X  The  following  LOOP  function  is  used  to  calculate  the  cdf  from  the 
x  negative  binomial  pdf  recursively. 


LOOP(L,  PTMP(I,J,K,L+1)$0KLVL2(I,J,K,L)  = 

Q(I,J)*C(NOBS(J>+ALPHA(J)+ORD(L)-l)/CORDCl)))*PTMP(I,J,K,L)  ); 


F/ 


& 


> 

o 

o 

H 

< 

Q 

ui 

u 

3 

3 

O 

E 

c. 

u 

E 


C(I,J.K)$0KLVL1(I, J.K)  = 

SUM(L$((0RD(L)-1)  LE  ( LMT< I , J )+ORD( K)-l ) ) ,  PTMPt I , J , K, L ) ) XPPC I , J ) ; 
FCI,J,K)$0KLVL1(I,  J,K)  =  LOG(  C(I,J,K)  ); 

DISPLAY  F; 

VARIABLES 

X(I,J,K)  1  if  stock  k  spares  for  component  j  of  end  item  i 
MOE  natural  log  of  probability  of  mission  completion; 

BINARY  VARIABLE  X; 

EQUATIONS 

WTCONST  weight  constraint  equation 

SELECT(I,J)  select  exactly  1  of  the  possible  stockage  levels 
MOEDEF  measure  of  effectiveness  definition; 

WTCONST..  SUM((I,J,K)$0KLVL1(I,J,K),  TWTC I , J , K)XXC I , J , K) )  =L=  WTCAP; 

SELECTCI, J)$(  0K( I , J )  EQ  1  ) . . 

SUMC  K$0KLVL1C I, J,K),  XCI,J,K)  )  =E=  1; 

MOEDEF..  SUM((I,J,K)$0KLVL1CI,J,K),  F( I , J , K)XX< I , J , K)  )  =E=M0E; 

MODEL  STOCK  / ALL/; 

OPTIONS  OPTCR  =  0.0015,  OPTCA  =  0; 

SOLVE  STOCK  USING  MIP  MAXIMIZING  MOE; 

PSUC  =  EXP(  MOE . L  ) 

DISPLAY  X. L, MOE. L, PSUC; 
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